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Abstract 



o : 

. A variety of coherent states of the harmonic oscillator is considered. It is formed by a par- 

^-H ' ticular superposition of canonical coherent states. In the simplest case, these superpositions 

> : 

■ are eigenfunctions of the annihilation operator A = Pid/dx + x)/^/2, where P is the parity 

m ■ 

' operator. Such A arises naturally in the q ^ —1 limit for a symmetry operator of a specific 

■ self-similar potential obeying the (7-Weyl algebra, AA'^ — q'^A'^A = 1. Coherent states for this 
I and other reflectionless potentials whose discrete spectra consist of N geometric series are 

^ • analyzed. In the harmonic oscillator limit the surviving part of these states takes the form of 

! orthonormal superpositions of canonical coherent states \e^a), k = 0,1, . . . , N — 1, where e 

^ ' is a primitive A^th root of unity, = 1. A class of g-coherent states related to the bilateral 

• • , g-hypergeometric series and Ramanujan type integrals is described. It includes a curious 

> ■ 

' set of coherent states of the free nonrelativistic particle which is interpreted as a q-algebraic 

. system without discrete spectrum. A special degenerate form of the symmetry algebras of 

self-similar potentials is found to provide a natural g-analog of the Floquet theory. Some 
properties of the factorization method, which is used throughout the paper, are discussed 
from the differential Galois theory point of view. 
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1. Introduction 



Replacement of commuting coordinate and momentum variables of a classical point particle by 
the operators x and p satisfying the Heisenberg commutation relation, 

[x,p] = xp — px = ih, (1.1) 

endows this particle with wave characteristics. According to the original definition, coherent 
states are the states in which corpuscular properties of a quantum particle are seen best of all. 
For the harmonic oscillator, such a qualitative motivation happens to be supported by the rich 
group-theoretical content built in the structure of coherent states, so that the symmetry ap- 
proach provides their alternative description. As a result, there have appeared several different 
quantitative definitions of coherent states which are equivalent only for the harmonic oscillator 
case. Independently of this non-uniqueness of definition, from both the physical and mathemat- 
ical points of view, coherent states of quantum mechanics are fascinating objects having useful 
applications in many fields |lll-@]. 

In the present paper we discuss coherent states associated with a specific class of one-dimensional 
Schrodinger operator potentials found in 0-0. The general class of these 5e//-szmz/ar potentials 
is defined with the help of g-periodic closure @, |^ of the dressing chain, or the chain of Darboux 
transformations. The latter transformations are known to be closely related to the factorization 
method Symmetries of the self-similar potentials are described by some polynomial oper- 

ator algebras of order (A^ is the period of closure) which play the role of spectrum generating 
algebras. For A^ = 1, 2 these algebras coincide with known g-analogs of the bosonic oscillator 
and su{l, 1) algebras 0. The q-coherent states, defined as eigenfunctions of symmetry operators 
which lower the energy, have many interesting properties. In particular, the algebras depend on 
the parameter g^, so that g = ±1 coherent states seem to be equivalent; but this is not so. For 
the g-Weyl algebra system, the limit g ^ — 1 exists only when the potential is symmetric, and 
then the corresponding coherent states are described by a particular superposition of canonical 
coherent states. This superposition has a universal form, it represents an example of the Titulaer- 
Glauber coherent states [|1^] which was constructed in |Tl[| from a different idea. Its multimode 



oscillator analog defines the particular entangled coherent states which have a phase difference 
equal to 7r/2. The name parity coherent states is suggested for these and other more general 



two-term superpositions of coherent states for which the parity operator plays a crucial role in the 
definition. 

For the limiting values of parameters corresponding to the harmonic oscillator potential, the 
raising and lowering operators of the self-similar potentials' symmetry algebras become equal to 
powers of bosonic creation and annihilation operators. In this limit, part of the g-coherent states 
degenerate into orthonormal superpositions of canonical coherent states |e'^a), k = 0,1, . . . , N — 
1, where e is a primitive A^th root of unity, = 1. These are natural generalizations of even 
and odd coherent states \T2\. Note that there are g-coherent states which do not survive in this 
limit. In another degenerate limit, when the potential vanishes but q remains arbitrary, one gets 
a unique set of non-trivial coherent states of the free particle system. 

It is interesting that symmetry operators of the self-similar potentials are defined with the help 
of a scaling operator - a special form of the squeezing operator. Due to this fact, wave functions 
of these systems resemble wavelets. From a degenerate form of the symmetry algebra, when pure 
dilatation by a fixed number q becomes physical symmetry of the Schrodinger equation, one finds 
a natural g-analog of the Floquet theory. 

Despite their simplicity, group-theoretical roots of superpositions of canonical coherent states 
comprise an intrinsic possibility for construction of complicated physical systems whose coherent 
states share common properties with these superpositions. The author considers the description 
of the relationship between self-similar potentials and superpositions of coherent states as an 
important physical result of this paper. However, a wider aim of the work is to discuss exactly 
solvable potentials and their coherent states from the general functional- analytic point of view on 
the basis of old and new examples. 

The paper is organized as follows. Before moving to the analysis of complicated situations, at 
the end of this section we give a brief account of the canonical coherent states of the harmonic 
oscillator. In section 2 we construct non-standard coherent states of the same system and consider 
their relation to the Titulaer-Glauber coherent states. In section 3 we discuss a universality of 
the derived superpositions of coherent states. An interesting set of coherent states of the free 
particle determined by the pantograph equation and its generalizations is considered in section 
4. The simplest potentials with g-deformed symmetry algebras and non-trivial discrete spectra 
are described in section 5, where some properties of the associated coherent states are analyzed. 
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In the next, 6th section we present a general hierarchy of Schrodinger operators whose discrete 
spectra consist of geometric series generated by the specific polynomial quantum algebras. Two 
particular systems arising from N = 2, g = — 1 and = 3, q = 1 closures of the dressing chain 
are considered in section 7. In section 8, a g- analog of the Floquet theory is outlined. Section 9 
contains a discussion of integrable potentials and coherent states from the differential ( "quantum" ) 
Galois theory point of view. Some concluding remarks are given in section 10. The paper has a 
formal character; we consider mostly theoretical aspects of the chosen (stationary) systems rather 
than their possible experimental implementations. Present analysis of coherent states for the self- 
similar potentials arose from the investigation of g-oscillator algebra at roots of unity performed 
T3| . The results of this work have been reported by the author in H^, |T^, and a part of them 



m 



has been published in p!5[| . 

In practical applications it is convenient to use the coordinate representation of ( |1.1| ), where 
p = —ihd/dx. For simplicity we use the system of units where Plank's constant h is equal to 1. 
In terms of the ladder operators a\ a: 

relation ( |1 . 1]) takes the form of the bosonic oscillator, or Weyl algebra: 

[a,a^] = l. (1.3) 

The number operator = a'^a satisfies the relations [A^, a^] = a\ [N,a] = —a. In appropriate 
units the Hamiltonian of a harmonic oscillator is equal to A^ up to a constant term: 2H = {a"!", a} = 
—(P/dx^ + x^. The energy spectrum and orthonormal eigenstates are: 

H\n) = En\n), E„ = n + 1/2, \n) = ^{a^Y\Q) ^ 

yn\ 

where the vacuum state |0) is defined from the equation a|0) = 0, (x|0) = vr^^/'* exp(— x^/2) (we 
set the phase of this state equal to zero). 

Coherent states of the harmonic oscillator may be defined either as eigenstates of the annihi- 
lation operator a: 

a\a.) = a\a), (1.4) 



or as a result of the application of the displacement operator to the vacuum: 

^ e^'^'-^'^lO) = e-l"l'/2 ^ ^\n). (1.5) 



=0 vn! 



Both definitions are essentially equivalent and give: 



V'a(x) = {x\a) = vr ^/^exp l^- — - (-^ - afj . (1.6) 

The states ( p. .41) are defined up to a phase factor expix(Q;, a*), where x{(^) ol*) is an arbitrary real 
function such that x(0,0) = 0. Only a special choice of x corresponds to ( |1.6| ). Note that the 
shift of X by a real constant Xq is a canonical transformation which is not completely equivalent 
to the shift of a: 

Since the bound states wave functions are expressed through the Hermite polynomials Hn{x), 

(:,\n) = ^(£Le-^'/^ (1.7) 



'2"n!A/vr 

the relations (|1.5| ), ( |1.6|) lead to the generating function for these polynomials 



n=0 



V2- 



The strongest quantitative measure of differences in the behavior of quantum and classical 
particles is expressed by the Schrodinger- Robertson uncertainty principle [|l^] : 



A = a^^app - alp > ^, (1.8) 

where abc = ^{bc + cb) — {b){c), and the angular brackets denote averaging over an arbitrary 
normalizable state for which the mean values are well defined, (6) = {ip\b\ilj)- 

Averaging over the coherent states \a) one finds: (x^) = | + (x)^, (p^) = | + (p)^, {px + xp) = 
2{x){p), where 

(x) = {a + a*)/V2, (p) = {a- a*)/tV2. 

In this case A = 1/4. Since axp = 0, this is the lower bound of the Heisenberg's uncertainty relation 
as well: (T^xCpp = 1/4. Note that the latter equality does not determine uniquely coherent states. 
One may scale x in ( |1.6| ) by a real number and preserve minimality of the product cr^^Cpp- The 
resulting states are called squeezed states; a more wide class of states corresponds to the lower 
bound of ( 11.81) . It is necessary to impose additional constraints in order to get (|L^) uniquely 
161, O. 
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The elementary example considered here illustrates three possible definitions of coherent states 
for an arbitrary system: 1) as eigenfunctions of some symmetry operators lowering the energy, 2) 
as an orbit of states generated by a chosen group element from a fixed state, and 3) as minimum 
uncertainty states for some physically significant operators. It is the first definition that we employ 
in this paper. 



2. Parity invariance and superpositions of coherent states 

The formulas (|1.2|) - ([L6|) are well known and widely used in quantum physics. First, we would 
like to find their analogs for a non-standard realization of the bosonic oscillator algebra. It is 
easy to see that this algebra has a nontrivial automorphism (i.e. a map onto itself), or canonical 
transformation associated with the parity operator P: 

PxP = -X, PpP = -p, P^ = 1, P^ = P. (2.1) 

Moreover, the transformation of x and p to the hermitian variables 

X = —ipP, p = ixP, 

is also canonical: [x, p] = i. Although this is a quite simple fact, it leads to the non-trivial reshaping 
of coherent states. Let us define new creation and annihilation operators: 

so that A + A' = \^x, A — = i\/2p. Evidently, the algebra, Hamiltonian, and vacuum state 
of the harmonic oscillator defined by the operator A coincide with the ones considered in the 
preceeding section. However, there is an essential difference in the structure of energy eigenstates 
generated by A'^: 

|n)„e«, = -^{A^r\0) = (-l)^ln), (2.3) 



[2A) 



where the sign factor (—1)'^" has the form: 

r 1, n = Ak, Ak + 1, = 0, 1,2, . . . 

I -1, n = Ak + 2, Ak + 3. 

This follows from the relation {A^)'^ = — (a"!")^ and the parity invariance of the vacuum. 



Denote by \o.)p the eigenstates of A, A\a)p = a\a)p, or 

{d/dx + x)il)^{x) = V2ailj^{-x), i^ai^) = {x\<^)p, (2-5) 

This is not an ordinary differential equation, but it can be easily solved using the relation A"^ = 
—a? . Picking out the appropriate combination of two linearly independent eigenstates of with 
the eigenvalue — a^, we find: 

\a)p = (e-'^/^\ia) + e'^'^\ - la)) , (2.6) 

p{P\a)p = {f3\a) = exp{f3*a -^\a\'- 

where |a) are the canonical coherent states. In the coordinate representation one has 

p \pi ( c? — \ol\^ — x^\ , r \ 
^aW = ^e^P| 2 I cos(V2q;x- -). (2.7) 

We call ( |2.7| ) the parity coherent states, because the parity symmetry plays a central role in their 
definition. As it is argued below, analogous states can be constructed for an arbitrary symmetric 
potential. Since \<y)p is defined by ( p..5| ) with A, A^ instead of a, , one can derive the following 
generating relation for the Hermite polynomials 

V ^ \ H^{x) = V2e'" cos ( 2tx - - ) . 

The |a)p's are not minimal uncertainty states for the variables x and p when a 7^ 0: 



a.. = (!-(«- a* f - (a + a*) V^I"l')/2, (2. 
= (1 + (a + a*f + (a - a*)2e-^l"l')/2, 
cr.p = («*'-«2)(l + e-^H^)/2^, 

so that 

A = l(l + p(l-(l + p)e-'')), p = 4|ap, 

and the minimum is reached only for p = 0. However by construction itself, the parity coherent 
states minimize the product of uncertainties in the new canonical variables x and p. 
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The minimum of a^x (|2.8| ) is reached for a = a* = 1/2: Cxxlmin — ~ ~ 0.32, i.e. 

there is a squeezing for small However, the states \<y)p differ from the squeezed states |T7 
arising as eigenstates of the annihilation operator after the canonical transformation 

A = S'^aS = cosh b| a + -- sinh \z\ a\ 



\z\ 



generated by the unitary operator 



S{z) = e(-'^^'-^*'^')/2. (2.9) 



If one multiplies this A by the parity operator from the left, then eigenf unctions of the resulting 
operator will be given again by superposition ( ^.6| ), but now with the |q;)'s on the r.h.s. being 
replaced by squeezed states. Note that in the Bargmann-Fock representation of the harmonic 
oscillator algebra, when A = Pd/dz, A^ = zP, one has V'q (-z) oc cos(q;2; — 7r/4). 

It is easy to construct analogs of ( |2.6| ) for the two oscillator algebra: [aj,(3.fc] = Sjk, [aj,ak] = 
0, j,k = 1,2. Again, the combinations Aj = Paj, where P is an operator which inverts both 
space coordinate axes, satisfy the same algebra. Such a transformation affects only the sign of 
the energy eigenfunctions. Eigenstates of the Aj-operators, Aj\ai,a2)p = aj\ai,a2)p, are given 
by the superposition: 

\ai,a2)p = (e~'''^'^\iai)i\ia2)2 + e'"/^] - iai)i\ - 102)2) , (2.10) 

where \a)j are the canonical coherent states of the j-th degree of freedom. This formula is obtained 
by choosing an appropriate linear combinations of the a^-operator eigenstates. Generalization of 
( |2.10| ) to an arbitrary number of oscillators is obvious. 



The unitary operator U that transforms ( |1.2D into ( ^.21 ) (or its multidimensional analog) is 



easily found due to the relation between the operators P and in the Hilbert space: P = 
(-1)^ = expivrA^. It is 

This is a special case of the canonical transformations generated by the operator 

U = e-'^^^\ U^aU = e'^^^^a, ip{N) = e{N) - e{N + 1), (2.11) 
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where ^(A^) is a non-singular function (such transformations were discussed in [jT8|, |T9| without 
relating them to superpositions of coherent states). In the classical case, when Hd = a*a, a = 
(x + ip) I -\/2, relations ( [2.11| ) correspond to a rotation of a through an angle depending on energy: 
Wall = e*'^("*")a which is obviously a symmetry of the system. 

Applying the operator W ( p. Ill ) to \a), one finds eigenstates of the operator A = e^'^^^^a: 



(2.12) 

n=o Vn\ 

These are the generalized coherent states of Titulaer and Glauber [l^. The states ( p.6|) belong to 
their particular subclass characterized by the periodicity condition 



9{n + M) = 9{n), 



(2.13) 



imposed upon the function 6{n) (in our example M = 4, cf. (|2.4|)). In this case, one has pO|, pT 



\coh) = ^ Ck\e^a), e = e 



2-Ki/M 



k=0 



i.e. \coh) is a superposition of M coherent states with the parametrizing variable a modulated by 
the powers of a primitive Mth root of unity. The states \e^a) are linearly independent, so that 
they can be orthonormalized. For the simplest M = 2 case these orthonormal superpositions take 
the form of even and odd coherent states fT^ 

\a) ± \ — a) 
'"^^ " V2±2e~2H^' 

The general root-of-unity analogs of these states \ai), / = 0, 1, 

following form 

M-l 



M - 1, {ai\am) = Sim, have the 



(2.14) 



m=0 



Mk+l 



Ma(a)e-i"i y - , 



Mk + l), 



\C,{a) 



/M-l 



-1 



M 



y e-'"^exp(e'"|an 



\m=0 



It is easy to derive these formulas using the kernel of the finite-dimensional Fourier transformation 



1 



Im 



dm 



M 



F^F = FF^ = 1, 
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F — F, {F )im — SioSmO + SM-l,ni, F — 1. 



The states ( p.l4|) do not belong to the class ( p.l2| ) because for arbitrary values of a they cover 



only a part of the Hilbert space, 

(fa \ai){ai\ = vr^, (fa = d{Re a)d{Im a)/iT, 



M-l 

whereas the Titulaer-Glauber coherent states form an overcomplete set of states. The possibility 
to split Hilbert space of the harmonic oscillator to an arbitrary number of orthogonal subspaces 
is related to the fact that the projection operators tt; are conserved, 

^'^ m=0 

The distinguished property of the parity operator is that it is the only symmetry of separate 
kinetic and potential terms of the harmonic oscillator Hamiltonian. Moreover, the definition 
Pf{x) = f{—x) works for any function independently of its membership in the Hilbert space (for 
complex X this is just the rotation of the complex plane by vr). It is only for the functions which 
can be expanded over the normalizable states \n) that the operator expinN coincides with parity, 
e.g. for the non-normalizable eigenfunctions of Hamiltonian the action of this operator is formally 
equivalent to multiplication by some phase factor which is not related to the transformation 
X —>■ —X. Therefore it is not clear what are the analogs of the operators ex])i6{N) beyond the 
Hilbert space context. 

Although the finite-dimensional truncation of the Titulaer-Glauber coherent states was dis- 
covered a long time ago |2y, |21]], only recently in |Tl| have Yurke and Stoler derived explicitly 



the superposition ( p.6|) , with ia being replaced by a. In their approach it emerged as a result of 
a fixed period time evolution of the standard coherent states governed by a specific Hamiltonian 
oc A^'^, k even. The squeezing and other properties of these and more general finite-term super- 
positions of coherent states have been analyzed in [E^-Pl]. The states (|2.1(J|) belong to the class 



of so-called two-particle entangled coherent states p2| , |25| , p6[] , which have the intrinsic property 
of non-factorizability into the product of one-particle states. For the construction of multimode 
analogs of the even and odd coherent states, see |^ . Recent interest in superpositions of macro- 



scopically distinguishable quantum states (Schrodinger cat states) like (p.6|) or the entangled states 
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like ( |2.10|) is inspired by new possibilities to create them with the help of optical techniques. This 
in turn provides interesting experimental tests of the basic principles of quantum mechanics (see, 
e.g. O, EH EBl EH] and references therein). 



The general set of coherent states constructed with the help of canonical transformations based 
upon the parity operator has the following form. Consider the unitary operator V, 

V = cosip + tPsmip, V^V = VV^ = 1, (2.15) 

where (f is an arbitrary parameter, and construct the Weyl algebra generators analogous to ( p.2|) : 

A = Va = aV\ = a)V^ = Va\ (2.16) 

[A, A^] = [a, a^] = 1, A^ = a^, {X^f = {a^f. 

The eigenfunctions of the operator A, or parity coherent states, are easily found 

= ^ ((1 + e-^^)^„(a;) + (1 - e-^^)^_„(x)) , (2.17) 

with ipaix) defined in ( |1.6|) . Calculating the uncertainties of x and p in these states, we obtain 



A(y.) = ^ (l + psin^ ^(1 - (1 + p)e~P)) , p = 4|a|2. 

It is seen that the choice Lp = tt/2, which corresponds to the Yurke-Stoler coherent states, is 
extremal - for it the value of A is maximally deviated from the standard coherent states case 
(f = 0, when A = 1/4. 

Let us discuss the Titulaer- Glauber states when the phases 9{n) satisfy the following q- 
periodicity condition: 

e{n + M) = qe{n), (2.18) 
which is a simple g-deformation of the condition ( p.l3| ). When q ^ 1 one can renormalize 6, 



6 = 6 — (f)/{l — q), and get 6{n + M) = 9{n) + 0, which is a quasiperiodicity condition. However, 
the effect of such a shift by is equivalent to the multiplication of a by the factor exp i(f)/M which 
is harmless for the representation of \coh) as a finite-term superposition of canonical coherent 
states. 

For generic values of q there is no split of \a) onto a superposition of a finite number of coherent 
states. For M = 1 one has 6{n) = and 

\coh) = e-HV2^ « = (2.19) 
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However, when the parameter g is a primitive Mth root of unity, = 1, the sum ( p.l9| ) is 
truncated: 

M-l 1 M-1 

\coh) = Y: Bi{<P)Wa), 5,(0) = V7 E g-'™e^^'"^. 

1=0 m=0 

The coefficients Bi are defined by the finite sums of exponentials which were already encountered 
in the calculation of normalization constants for the orthonormal states \ai) (|2.14|) . We would 
like to note that the g-deformation (|2.18| ) is not related to the g-coherent states to be discussed 
below, although the discrete energy spectrum of the corresponding Hamiltonians is found from 
the similar formula E^+m = Q'^En- 



3. Universality of superpositions of coherent states 

Let us discuss whether the superpositions of a finite number of coherent states considered above 
are characteristic only for the harmonic oscillator case or their form carries universal character 
applicable to any system. Unfortunately, there is no completely satisfactory definition of coherent 
states la) for an arbitrary Hamiltonian 

H = -^^+u{x), (3.1) 

even if the potential u{x) is an analytical or infinitely differentiable function of x (i.e. when the 
random, singular, and potentials from C^, k < oo, are excluded). Denote by |A) the physical 
eigenfunctions of the Hamiltonian: 

H\X) = E{X)\X), 

where A is some index labeling the spectrum; A = n = 0, 1, . . . for the discrete eigenvalues (the 
ordering E{n) < E{n + 1) is assumed), and for the continuous spectrum A may be thought of 
as some continuous variable such that the spectrum -E'(A) is monotonically covered by variation 
of A, say from Aq to oo. We assume that the continuous spectrum states are normalized by the 
condition (A la) oc 6{X — cr). From the completeness of generalized eigenfunctions of H in Hilbert 
space, 

Y^\n){n\+ dX\X){X\ = l, (3.2) 

n=0 -^^0 
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where rih is the number of bound states, it follows that for any definition of coherent states 
these states can be expanded over |n) and |A): 

" i) ~ ^ POO 

\a) = ^ Cn{a) \n) + dX cx{a) |A). (3.3) 

n=0 •^'^0 

The continuous spectrum exists whenever, for x oo (or — oo), the potential is not bounded 
from below or it is bounded from above. In these cases one cannot exclude in general the second 
term in expansion ( ^.3] ). Such possibilities are rarely discussed in the literature; some models of 
coherent states built only from the continuous spectrum states are described in the next three 
sections. 

Suppose for a moment that the Hamiltonian H has only a discrete spectrum, i.e. rif, = Aq = oo. 
Then, following many existing examples it is natural to assume that coherent states \a) 
represent a generating function for stationary energy states of the form 

oo 

\a) = J2a''cn{\a\)\n). (3.4) 

n=0 

Since the coefficients c„ depend only on the modulus of a, these states are complete, 

/oo « 
(i^a;p(|a|)|a;)(a;| = ^ |?T,)(m| / cfap{\a\)a*"^a^CnC^ = 1, 

n,m=0 

provided there exists a measure density p(|a|) which satisfies the relations 

/•oo 1 

/ rfrp(r)r2"+i|c„(r)|2 = - 
Jo 2 

for arbitrary n = 0, 1, . . .. When the variables of c„(r) separate, c„(r) = Cng{r), which is typical 
for the ladder operator approach, this is a moment problem. Theoretically it is possible that p(r) 
is defined non-uniquely in which case there are many physically distinguishable representations of 
observables in the coherent states basis. 

An interesting definition of coherent states based on the uncertainty principle was suggested in 
for a wide class of potentials. It uses the fact that for a classical particle moving in a convex 
potential one can make a non-canonical non-linear change of phase space variables such that in 
the new "coordinates" the particle's dymanics is described by the harmonic oscillator equations of 
motion. After quantization, the minimum uncertainty states of these harmonic motion "position" 
and "momentum" operators are called coherent states. By definition, this procedure is tied to 
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the quasiclassical approximation. Its general group-theoretical meaning is not clear to the author, 



probably the approximate dynamical symmetry approach of [IIH] can be useful in this context. 



Another constructive definition of coherent states for generic discrete spectrum systems has been 



suggested by Klauder In this approach, the requirement that time evolution of coherent 

states be equivalent to the change a e*'^*a is taken as a basic property. Such states do not 
spread but for them one has more complicated expansions than (|3.4|) . 

In the following we utilize a different definition of coherent states; namely, we assume that 
they are defined as eigenstates of some lowering operator A, i.e. the operator which maps a part 
of physical solutions of a given stationary Schrodinger equation to the physical ones with lower 
energy. Even for systems with purely discrete spectrum this requirement does not necessarily 
mean that A annihilates the ground state or that A is the lowering operator mapping discrete 
spectrum eigenfunction \n) to the closest from below state \n — \) (it may jump over some physical 
states which is always so for the continuous spectrum). For example, similar to the situation in 
0, 1^ , y4 may play the role of both lowering and raising operators for different ranges of energy. 

We use the factorization method p|, Q as a tool for searching for such a symmetry operator A. 
Let us factorize the Hamiltonian (|3.1|) , i.e. represent it as a product of two first order differential 
operators conjugated formally to each other: 

H = a^a + EQ, a = d/dx + f{x), = —d/dx + f{x), (3.5) 

where Eq is some constant. The potential u{x) and superpotential f{x) are related by the Riccati 
equation, u{x) = P{x) — f'{x) + Eq. The zero mode of a, 

aV'o(x) = 0, Vo(x) oc e-r^(^)'^^ 

is the generalized eigenfunction of the Hamiltonian with the eigenvalue Eq. When this function is 
normalizable and nodeless, Eq is the ground state energy. If a is a symmetry operator, i.e. if it 
maps physical eigenstates of the Hamiltonian H onto themselves, then coherent states could be 
defined as eigenstates of a, but there may be non-uniqueness even in this simple situation. Indeed, 
the factorization involves an arbitrary unitary operator T: 

A^A = a^a, A^ = a^T, A = T'^a, = T-\ (3.6) 

Actually, one may write A^ = a)C, A = Da, where C and D are two operators satisfying CD = 1. 
The operators A^ and A are conjugated to each other if C = D'', i.e. when D is an isometric 
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operator Only the additional requirement DD^ = 1 makes D unitary, we restrict ourselves 



to this case. 

Coherent states of potentials for which A = T~^a is a symmetry operator, for some unitary 
operator T, are thus defined by the equation: 

^o(a;)-^74~V = OiTijj^ix), (3.7) 
ax ipo{x) 

where ipo{x) is an eigenstate of the Hamiltonian (not necessarily a physical one). The natural 
extension of this definition involves on the l.h.s. of (|3.7|) a differential operator of the A^th order. 
In that case A^A is equal to an order polynomial of a Hamiltonian, i.e. one has a generalized 
factorization scheme which will be described in section 6. 

Suppose now that the expansion ( p.4|) takes place (it does not mean that there is no continuous 
spectrum, there may be accumulation points such that the set \n) is closed under the action of A). 
Then, one can introduce the formal number operator, A^, satisfying N\n) = n\n). Acting upon 
such |a) by the unitary operator W = e^^^^\ one gets the Titulaer-Glauber type coherent states 
for the chosen class of potentials. The periodicity condition 9{n + M) = 9{n) leads again to the 
finite term superpositions of \e^Q). Therefore the form (but not the normalization constants) of 
the superpositions \ak) (|2.14|) is universal - for any H they perform a split of discrete spectrum 
Hilbert subspace onto orthogonal components. Note however, that the symmetry properties of 
these superpositions may be different; in particular, the general even and odd coherent states are 
not eigenstates of the parity operator. For instance, for the shifted harmonic oscillator potential 
u{x) = (x — XoY one has: 

(x|a±) oc e-("-"°)'/2 |^gV^{x-xo)" ± g-V2(x-xo)a^ ^ ^3^8) 

which are eigenstates of the hermitian conserved charge Pexp(2xoci/(ix). For the general asym- 
metric potential, are eigenstates of the abstract operator expivrA^, which does not have a 
simple form in the coordinate representation. 

Similarly, one can always define upon the discrete spectrum an analog of the Yurke-Stoler 
coherent states, but they will be eigenstates of the operator Pa, an analog of ( |2.2| ), only in the 
case when Pa = —aP, which assumes that the potentials are symmetric, u{—x) = u{x) (if one has 
Pa = aP, then the symmetric and antisymmetric eigenfunctions of a diagonalize simultaneously 
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the operator Pa). We come thus to the conclusion that the notion of "parity coherent states" is less 
universal than that of even and odd coherent states and their higher root-of-unity generalizations 
defined by the abstract formula ( |2.14 ). The special name for the states ( 2.17| ), and similar ones, 
was introduced in order to distinguish Titulaer-Glauber states with the M = 2 periodic phases for 
symmetric and asymmetric potentials (in the first case there is no need for the abstract operator 
N). 

Single-valuedness of the expansion (|3.4| ) under the rotation a e^'^^a is the key property 
allowing to build finite-term superpositions of coherent states ( 2.14 ). In general such property 
does not hold. In the next section we consider a model where only the continuous spectrum piece 
is present in (|3.3| ) with c\{oi) oc a'^"^^!^^^^ . In this case the change a e^a, e'^^ = 1, does not 
provide a split of Hilbert space onto the finite number of orthogonal components. 



4. Coherent states of the free particle 

The above definition of coherent states ( |3.7| ) is applicable to systems with continuous energy spec- 
tra. Consider the simplest possible case of zero potential, H = —(P/dx'^, for which the generalized 
eigenfunction of lowest energy is a simple constant. Solutions of the Heisenberg equations of 
motion x = pot + xq, p = Po, where t is the time variable and ^052:0 are operators at t = 0, are 
identical with the classical ones. If this coincidence would be taken as the basis for the definition of 
coherent states, then any normalizable state of the free particle has to be considered as coherent. 
However despite the same form for the equations of motion, even for the (po) = case, when 
the classical particle stays at the point xq, the quantum particle tries to occupy the whole space: 
Cxx oc (Po)^^) t ~^ s-^d non-spreading wave packets do not exist. 

The simplest factorization of the free particle Hamiltonian is obvious and we have a = d/dx. 
Let us find eigenfunctions oi A = T'^d/dx when T is a unitary operator performing an afiine 
transformation, or the parity operator. If T = 1, then ipa{x) oc e"^ which are bounded (but 
unnormalizable) functions for purely imaginary a: a = ip, —00 < p < 00, so that ipip{x) coincide 
with the momentum eigenstates. If T is a translation by h operator, then 

dipa (x) /dx = aipa + h) 

and again there are simple solutions of the form ipa{x) oc e*^^ but now a has a real part: a = 
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ip^-wh^ If T = P, the parity operator, then one gets superposition (|2.6|), ipl^{x) oc cos(ax — vr/4), 
a real. Since these states are not physically reahzable (they belong to the continuous spectrum), 
their "coherence" is formal. This is related to the fact that the chosen symmetry operators A 
are integrals of motion commuting with the Hamiltonian. As is seen from the considerations 
given below, such A's appear from real lowering operators in a special limit such that, in fact, all 
Hamiltonian eigenstates acquire a flavor of coherent states (this is reminiscent to the approach of 
1^, where coherent states are defined as eigenstates of some combinations of integrals of motion). 
Examples of symmetry algebras for which the Hamiltonian eigenstates might be counted among the 
coherent states are given below. In the corresponding cases the operator A may be simultaneously 
the lowering and raising operator, and an integral of motion for different parts of the spectrum. 
Let T be the scaling operator: 

Tfix) = ^/Wlfiqx), (4.1) 

where q is some real parameter, < < 1. For positive q, T is just the squeezing operator in a 
special form: T = S{z = Ing) = where a^,a are given by ( |1.2|) ; and for negative q it 

is a product of the same S and the parity operator: T = SP. For simplicity we assume in this 
section that q is positive. Now the operator A = T'^a is not an integral of motion, but the raising 
operator for the A > solutions of the Schrodinger equation 

H%l){x) = -iIj"{x) = Xipix). (4.2) 

The commutation relations 

AA^ = q^A^A, AH = q^HA, HA^ = q^A^H 



look similar to some of the defining relations of quantum groups [p2| , [33|| . In fact, they are the 
progenitors of g-deformed oscillator algebra (see below). Eigenfunctions of the operator A are 
determined by the differential equation with deviating argument 

tjj'^{x) = ay/qipa{qx), < g < 1. (4.3) 



The initial value problem for ( [4.3| ) is qualitatively different from that for the ordinary differential 



equations because the initial conditions have to be fixed now on the interval [gxo,Xo]. When xq 
is a fixed point of the scaling transformation, i.e. = or oo, this interval shrinks to one point. 

18 



When xq = oo, one fixes solutions by taking the asymptotic form of ipa{x) from some class of 
permitted functions |0. For xq = it is natural to impose the initial condition ipaiO) = 7 < oo. 
Then equation ( ^.3| ) has the unique analytical solution: 

00 n(n—l)/2 

M^) = 7 E - — —i^^V^^r, (4.4) 

n=0 ^• 

which is an entire function of x for any finite |a|. However as is shown in (see also |]35[), the 
X ^ 00 asymptotics of any solution of ([4 .31) is dominated by the factor exp(— In^ x/ In g^), i.e. all 
solutions grow at infinity so that the functions ipa{x) do not describe physical states. Actually, 
this could be expected from the fact that we were diagonalizing the raising operator. 
Consider eigenfunctions of the lowering operator , 

ax 

or 

ij'^ix) = - aq-^/^ij^{q-\), < q < 1. (4.5) 

Formal series solution of this equation with the boundary condition ip{0) = 7 looks similar to ( [4. 41 ) 
with q being replaced by q~^. But the radius of convergence of this series is equal to zero, i.e. 
there are no solutions analytical at zero. This does not mean that there are no solutions at all, 
in it was shown that in fact there are infinitely many non-analytical solutions of ( [4.5|) from 
C°° satisfying boundary condition tp{0) = 7. Moreover, for any a there are solutions with the 
|x| 00 asymptotics oc exp(ln^ |x|/lng^), i.e there are functions "ipaix) which are normalizable. 
One can expand such tpa{x) over the basis of Hamiltonian eigenfunctions 

Mx)= / dpe'P^Mp), (4.6) 
Jo 

which is the positive momentum part of the standard Fourier integral. Similarly one can consider 
the expansion over e~*^^ (for negative q the two regions of p should be considered simultaneously). 
Substituting this expression into ( |4.5|) and solving the corresponding finite-difference equation for 
the formfactor (f)a{p), we find 

M^) = r dp e'^^hip) exp 1^^!^/!^, (4.7) 
Jo 2 in (7 
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where h{p) is an arbitrary function periodic on the logarithmic scale, h{qp) = h{p), normalized by 
the condition 

°° , , / m9 , ^ /■°° , / m9 In^p/lal - (argia)^ 



/o In q 

Due to the freedom in choice of h{p) these states can take various forms. 
Let us take, for example, the following h{p), 

oo 

hip) = a Y: q^'Sip-bq^), 

n=—oo 

which does not correspond to normalizable ipai^)- Despite the formality of the consideration in 



this case, we obtain the free particle's "coherent states" in the form of Dirichlet series (cf. ^ 



M^) = aexpi^^^ £ (-Y^'^"^, (4.8) 

21ng \ia) 

which defines a function bounded for any real x and < |a| < cxo. This is an example of the 
bounded but non- normalizable solution of ( [4.5|) . The derived expansion could be interpreted 
as a g- Fourier series for ipaix) because the sum goes over the trigonometric functions whose 
argument is modulated by the powers of q (after renormalizations, Cn exp(i5;g") may become 
the standard Fourier series in the limit g — > 1 due to the relation (g" — l)/(g — 1) ^ n, in our 
case this gives a divergent series). This situation does not seem to be related to the g- Fourier 
transformation considered in it resembles more the wavelet transform or the generalized 
Taylor expansions for atomic functions . 



In the above factorization of the free particle Hamiltonian we took £^0 = 0. A much more rich 
situation arises when Eq = — where /5 is a nonzero real number. This gives H = a)a — jS'^, 
a = d/dx+/3. Evidently a is an integral of motion, but from the properties of affine transformation 
operator T, 

Tf{x) = ^f{qx + l), 0<g<l, -cx)< / < oo, (4.9) 

where g and / are fixed parameters, it can be seen that A = T~^a is the raising operator for 
positive energy states. In order to save space we shall describe the whole hierarchy of such 
symmetry operators at once. Let us introduce the operators 

N / J \ N 



(4.10) 
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where (3^ are arbitrary real positive constants. It is easy to check that and A satisfy the 
following nonlinear algebraic relations 

N N 

A^A = n + Pt) , AA^ = I[ {<l'H + f3l) , (4.11) 

k=l k=l 

HA^ = q^A^H, AH = q^HA. (4.12) 
For = 1 these relations define a g-analog of the Weyl algebra, or g-oscillator algebra (see, e.g. 



AA^ -q^A^A = u}, u} = (3l{l-q^). (4.13) 



For N = 2 one gets a g-analog of the su{l, 1) algebra in the form considered, e.g. in p3| , p[ . For 
N > 2 one has polynomial quantum algebras |]^, . Note that in the limit g ^ 1 one does not get 
non-trivial algebras since the operators A and A^ start to commute, still being the differential- 
difference operators for / ^ 0. For g 7^ 1 it is possible to set / = by going to the reference frame 
where a fixed point of the affine transformation is taken as the zero point, we shall assume this 
choice below. 

Discrete series representations of the derived algebra are constructed by the action of the 
operators A and A'^ upon the Hamiltonian eigenstates with A 7^ 0. The lowest weight series have 
the form 

A^^^\x) = 0, V^f )(x) cx e-^-\ = 1, 2, . . . , iV, 

(x) cx (At)>f (x), H^i'\x) = -plq^-i;^-\x), n = 0, 1, . . . , 00. 

The wave functions 'ipk^\x) are not physical, their eigenvalues accumulate near the A = point 
from below. One can take the (3k in ( [4.10| ) to be purely imaginary numbers. Then 'ipl^\x) describe 
continuous spectrum states whose eigenvalues accumulate near the zero from above. In this case 
A^ is not a hermitian conjugate of A but still one can use representations of the algebra ( |4.11|) , 
( [4.12| ). The highest weight representations appear as follows (in the same notations): 

AVf (x) = 0, #(x)oce'^'='^-'^ 

^i")(a:) oc A>f (x), i/^f (x) = -^3lq~'--'^pl:'\x). 

The eigenvalues of these states are unbounded from below for real f3k and they go to infinity for 
imaginary f3k- In any case operator A is the lowering operator for negative A eigenfunctions, but 
it raises the energy of continuous spectrum states. 
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Physical eigenstates of the Hamiltonian H = —(P/dx'^ have the form 

\/ Ztt 



where a, a' = ±. The algebra generators act upon them in a simple way 

N 



Aij^{x) = q-'/' n {±^^ + Pk) (4.14) 



k=i 

N 



A^ij^ix) = n {T^qv^ + Pk) i^t^ix). (4.15) 



k=l 



We conclude that the free particle's Hilbert space provides a unitary realization of the quite 
complicated symmetry algebras. In the following sections we describe a generalization of this 
construction to nontrivial potentials along the lines of . 

Consider coherent states of the above algebras defined as eigenstates of symmetry operators 
lowering the energy. Let us analyze first eigenstates of the "annihilation" operator A, Aipai^) = 
aipa{x). For the g-oscillator algebra such a definition has been considered, e.g., in [^, 

|48| , ^ and many other recent papers. In our model these coherent states for = 1 are 
defined by the differential-delay equation 

?/'^(x) = ay/qi)a{qx) - (3ii)a{x), (4.16) 



known in the literature as the pantograph equation |35|. Note that the initial value problem 



for this equation is highly nontrivial. Using the results of the detailed analysis of ( |4.16| ) given in 



3^, it is possible to see that solutions having finite fixed values at x = are normalizable 

near the x = oo point when |a| < (3i but for x —oo they diverge exponentially fast. A similar 
situation holds for the general > 1 symmetry algebras. We conclude that there are no physically 
acceptable coherent states of this type in the free particle model. This is caused by the absence 
of a discrete spectrum. In the more complicated realizations of (|4.11|) , ([4.12|) such states do exist. 



A different type of coherent states for the g-Weyl algebra has been constructed in [|I5|. These 
coherent states are defined as eigenstates of the operator A\ which lowers the energy of the A > 
states, A^ipctix) = aipa{x)- For the free particle realization of the g-oscillator algebra these states 
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are determined again by the pantograph equation but, now the dilation parameter is bigger than 
1: 

iIj'^{x) = -aq-'/^Mq~'^) + /^ig^Va (4.17) 

This time there are infinitely many normalizable functions ipa{x) for |a| > decreasing as 
x'^, q"^ = a/q^^'^Pi, for |x| oo. All of them have finite values at x = but they are not 
analytical near this point (i.e. Taylor series expansion does not converge). For arbitrary N, 
coherent states of this type are defined by the generahzed pantograph equation |^ 



n ^ + f^kj Mq^) = aq-'/^ij^{x), (4.18) 

whose solutions are not expressible in terms of the classical special functions. Expansion of 
normalizable solutions in the basis of Hamiltonian eigenfunctions is given by the integral 

,(,) _ C(a) A^^e^^^^c/A 



<l(a;) = / — ^ , (4.19) 



Inp/ag^/^ + 27ris 
ds = : — 5 , p = Pi---Pn, 

where s = 0, ±1, ... is an integer enumerating linearly independent states and it is assumed that 
< arga < 2tt. In ( ^.19| ) and below we use the standard notations for g-products 

(a;g)oo = n(i -ag''), {a;q)n = ■ 

For positive integer n one has 

(a; q)n = n (1 - («; = \ , X • 

The normalization constant C{a) is given by the integral 

This constant is finite for \a\ > p, which is the region of definition of the coherent states. For any 
the states ( [4.19|) have the common qualitative feature of non-analyticity near the x = point. 
The origin of this property and possible physical consequences deserve further investigation. 
One can look also for solutions of (|4.18|) in the form of series similar to ( [4.8|) . For example, for 



= 1 one can write: 



^„(x)oc f: (f^) (^^g/A;g)„e'^^"^ (4.21) 



a 
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where 6 is an arbitrary constant, q < 6 < 1. These non-normahzable solutions are bounded for 
any x provided |«| > g^^^/Si, however it is not clear whether they should be taken into account in 
order for coherent states to be complete. 

It is worth mentioning that the pantograph equation appears in various problems. It has been 
encountered in a description of the light absorption by interstellar matter the collection of 
current by the pantograph of an electric locomotive |5^, some number theory problem, etc (for a 



list of applications, see ||3^). Here we have described another physical application of this equation 



as the one determining the free particle coherent states within the ladder-operator definition 
context. 

Two linearly independent A = solutions of ( |4.2|) can be represented in the form 

a ^ 1 

V'i(x) = l, ^p,(x) = x + -^J2ir- (4.22) 

^ - ? fe=l Pk 

Only the first one is bounded and belongs to the continuous spectrum of the free particle. Formally 
upon ipj{x), both operators A and A'' are diagonalized simultaneously: 

AiPjix) = pq^/'^-^ipjix), A^-^jix) = pq^-^/'^'^j{x). (4.23) 

These are the simplest models of the c-number "condensate" representations of quantum algebras 
discussed in [|13|, p3| . 

In a similar manner one can consider the free particle on the half-line. In order to make 
Hamiltonian self-adjoint it is necessary to impose boundary conditions at x = 0. The scaling 
operator defines physical symmetry only in the special cases ip{0) = or ip'{0) = Then one 



can define again coherent states as eigenstates of the lowering symmetry operators, but we shall 
not consider them here. 

A curious fact is that for = 1, q = i, Pi = one arrives at a simple differential- 

difference operator realization of the fermionic oscillator algebra. Indeed, it is not difficult to 
check that upon the states 

A\0) = 0, |0)oce-^/^, |1) = At|0) oce"*^/^ (4.24) 

the relations A^l) = and A\l) = |0) are satisfied, which means that AA^ + A^A = 1, A^ = 
(y4^)^ = 0. It is not clear whether this formal two-dimensional representation has a unitary setting. 
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The scaling operator T also provides an interesting possibility of a non-standard realization of 
the ordinary bosonic oscillator algebra. Let Tf{x, y) = f{qx, qy) and z = x + iy, then the relations 

[A, A^] = 1, A = T-^-^, A^ = zT, (4.25) 

define some deformation of the Bargmann-Fock realization (when q is complex the operator T 
scales and rotates z, we assume that < g < 1). The measure density p(|-2p) in the scalar 
product 



(^i|V^2) = J d'zpi\z\')Mz)Mz) (4.26) 
is found from the requirement for to be a conjugate of A. This gives the equation 

^^-,-V(l), (4.27) 

which we have just encountered in ( [4.5|) . Its solution satisfying p(0) = const can be represented 
in the form 

p{t) = / dp h{p)e~^^ exp — , h{q^p) = h{p). 

Jo 4ing 

Acting by A^ upon the zero mode of A, Ai/jq^z) = 0, ipo{z) = const, one finds the Fock space basis 

vectors: 

Mz) = Cn{q)z'' oc {zTri. (4.28) 
The normalization constants C„ have the form 



\Cn\'= ^, E (-l)%exp^ , (4.29) 



where hk are arbitrary constants appearing due to the nonuniqueness of the measure. More 
precisely, hk are the coefficients of the Fourier expansion of an arbitrary periodic function entering 
the measure h{p), 

Hp) = 1^ hkexp— . 

, In (7 

fe=— oo ^ 

The relation between hk and C„ shows that the moment problem for this measure (i.e. determina- 
tion of p{t) from given C„) does not have a unique solution. Physical applications of the described 
realization of the Heisenberg-Weyl algebra are not known. 
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5. Coherent states of the ^'-deformed harmonic oscillator 
potential 



In section 2 we gave a simple derivation of the Yurke-Stoler states ( p.6[ ) on the basis of a canonical 
transformation associated with the parity operator. Actually, it was inspired by the analysis of 
coherent states for the g-oscillator algebra, 

AA^-q^A^A = u, [A,uj] = [A\uj] = Q, (5.1) 

in the realization described in p . Let us consider this system in more detail. One can check that 
the pair of formal operators: 

A = T'\d/dx + f{x)) , A^ = {-d/dx + f{x))T, (5.2) 



where T is the scaling operator ( |4.1| ), satisfy ( |5.1| ) provided f{x) is a solution of the equation 

^ (fix) + qfiqx)) + fix) - q'fiqx) = uj, (5.3) 

derived in as the simplest self-similar reduction of the dressing chain for Schrodinger equation. 
The Hamiltonian of this system, 

H = A'' A -V = -d^/dx"^ + m(x), 

u{x) = fix) - fix) - z/, u = - q^), (5.4) 

satisfies the relations: 

AH = q^HA, HA^ = q^A^H. (5.5) 

For q = 1 one has f{x) = ujx/2, i.e. the standard harmonic oscillator. Suppose that f{—x) = 
—f{x), which corresponds to the symmetric potential. Then, in the limit uj —>■ the solution f{x) 
analytical at x = vanishes due to the initial condition /(O) = 0, and we get the zero potential 
model considered in the preceeding section with the factorization constant Eq = (the Eq ^ 
case corresponds to the symmetric solution f{x) = \fv). 

The analysis of [|I^ shows that for complex values of x and g, a solution of ( |5.3| ) analytical 



near x = exists and it is unique provided |g| < 1 or g is a primitive root of unity of odd degree, 
^2fe+i _ For (f'^ = 1 a solution may exist only for special initial conditions which, however, do 
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not guarantee uniqueness. If |g| = 1 but 7^ 1 then f{x) = ±yV are the only analytical solutions 
known to the author. For < |g| < 1 the function f{x) cannot be expressed in terms of known 
special functions. When g^'^+i = 1, and in the restricted case of g^'^ = 1, the problem is solved in 
terms of particular hyperelliptic functions characterized by the presence of additional symmetries 
of the lattice of periods. Since in the simplest cases, = 1, q"^ = 1, q is just the modular 
parameter of elliptic functions, the function f{x) (and its generalizations to be described below) 
comprises hidden "second" g-deformation properties of hyperelliptic, or finite-gap potentials. 

In the system (^)-( |0| ) has been derived from a special quantization of a simple model 
of classical mechanics (a particle in a finite depth potential) characterized by a quadratic Poisson 
algebra. In this picture equation (|5.3| ) has the form 

^-r- f fix) + e^'^f(e^'^x)) + f(x) - e^^V^le^^x) = c(l - e^^"), (5.6) 
ax ^ ' 

where Ti is Plank's constant, q = e^^, i] and c are parameters of the classical potential. Expanding 
this equation over h one finds successively the classical, quasi-classical, and so on, approximations 
to the exact solution f{x) (note that in this approach the potential is an infinite series over h, 
there is a large nonuniqueness, etc). 

Let X and be real, u > 0, and f{x) antisymmetric, f{—x) = —f{x). For g^ > 1 the function 
f{x) has singularities so that is not conjugated to A and the realization of ( p.l| ) is not 
unitary. For < g^ < 1 the function f{x) is bounded and has only one zero f^. These are the 
crucial properties sufficient for A'' and A to be well defined operators in the Hilbert space. In 
particular, the zero mode of A, A\0) = 0, (x|0) oc exp{— f{y)dy), is normalizable and describes 
the ground state of Hamiltonian H. As a result the whole spectrum of H is found from unitary 
representations of the algebra (|5.1| ). 

It is not difficult to see from ( [5. 5] ) that the spectrum may consist of three parts: a discrete 
one, describing bound states accumulating near the zero energy level from below, a continuous 
part going from zero to infinity and corresponding to scattering states, and, finally, a zero energy 
piece. The discrete spectrum is described by the lowest weight discrete series generated by 
from the vacuum state |0): 

|n) = ^^|0), A\0)=0, {n\m) = 6nm, 
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[n]\ = [n][n-l]\, [0]! = 1, [n] = {1 - q'^) / (1 - q') , 



1/2 



1 _ /72{n+l) , / 1 _ „2n 

|n) = E„\n), E„ = - i^q'". (5.7) 



It consists of one geometric series. Since zero modes of A are determined by the first order 
differential equation, it follows that ( |5.7| ) are the only physical states for E < 0. Indeed, suppose 
that we missed one physical state \E) with energy E < 0. Acting by powers of A upon \E) we 
get a sequence of states of lower energies. Since the potential is bounded this series should be 
truncated which is possible only if \E) is annihilated by some power of A, i.e. if the state \E) 
belongs to the series (p.7[). 

The same argument shows that the states with positive energy, E > 0, appear in the form of 
a geometric series infinite in both directions, 

H\n)x = Ag2"|n)A, n = 0, ±1, ±2, . . . (5.8) 
A\n)x = ^Ju + \q^^\n-l)x, 
A^\n)x = y^i/ + Ag2{n+i)|n + l)A, 

where A > is an arbitrarily chosen eigenvalue of H. (This representation of the g-Weyl algebra 
was discussed in |13|, cf. also [Q; note that it is not defined for g 1.) Theoretically 



these states could be normalizable, but in our case this is not so - they belong to the continuous 
spectrum. Indeed, from (|5.7|) , ( pl8|) it follows that the scaled potential q^u{qx) has the same 



spectrum as u{x) except of the lowest state with the energy Eq = —v. Similarly, q^^u{q^x) does 
not have the k lowest states. Taking the limit k —>■ oo, one gets a system without a negative 
spectrum, whereas the states (|5.8|) are not washed away. From the boundedness of the initial 
potential it follows that q^'^u{q^x) ^ for A; ^ oo because q^ < 1. This means that our potential 
is reflectionless, being obtained by a special infinite step dressing of zero potential. But for zero 
potential the series ( ^.8| ) with = (removal of the levels is equivalent to rescaling uj q^^uj 0) 
correspond to the continuous spectrum (see the preceeding section). A rigorous proof of the 
absence of positive energy bound states requires an estimate of the asymptotics of the potential 
57[| . The |x| ^ oo asymptotics of ( ^.4|) proposed in decreases sufficiently fast in order to 
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guarantee the absence of such exotic states, u{x) h{x)/x^ + 0(l/a:^), where h[qx) = h{x) is a 
bounded function. 

Upon the zero modes of Hamiltonian the operators A and behave hke integrals of motion, 
i.e. they commute with H. Since the relations A'' A = AJ^ = v can be satisfied by arbitrary 
invertable matrix, the dimension of this representation is not restricted. It may be either infinite- 



dimensional under additional requirements or just one-dimensional. In the latter case the 



creation and annihilation operators degenerate into complex numbers |T3|, ^ 



E\d) = 0, A\d) = v^e^^VO, A^\cl) = V^e'^\d), (5.9) 

where we assume that the state \d) is normalizable. Note that the limit g — > 1 is not defined. 
In our case these "classical" states correspond to the boundary between discrete and continuous 
spectrum. As we shall show below, in the g-oscillator model ( [5. 21 ) with /(O) = the corresponding 
wave functions are not bounded and their eigenvalues differ from those in (|5.9| ). 

Remark. Actually, the existence of the non-zero c-number representations is not very rare for 
quantum algebras. For example, for the Cartesian version of the slq{2) algebra 



describing dynamical symmetries of some discrete refiectionless potentials |59|, one can set J; 



l/(g ^ — q) and the algebra is satisfied. Such exotic representations are often skipped in discussions 
of applications of quantum algebras. 

As we have discussed already, coherent states of the first type of the algebra ( ^.1|) are defined 
as eigenfunctions of the annihilation operator A. These states are built from the lowest weight 
discrete series representation ( p.7|) (see, e.g. 



a" 

A\a,q) = a\a,q), \a, q) = C{a) \n), (5.10) 

n=o Juj"'[n]\ 

where 

\C{a)\ =e,2iz) = ^ = , z= , 

n=0 [q ^q )n 1^1 q )oo ^ 

is a g-analog of the exponential function, or i(y9o(0; g^, 2;) basic hypergeometric function. The 
general g-series of such type is defined as follows [Q: 

^ 62, &s ' ' / ^0 (g;g)n(&i;g)n- • • (&.;g)n 
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where r and s are arbitrary positive integers, and Oi, . . . , a,., fei, . . . , 6^ are free parameters. The 
states ( [5.10| ) are normahzable only for lap < u. Since the Hilbert space of our model is larger 



than the discrete spectrum Fock space spanned by \n), the states ( p.lUD are not complete. There 
remain two parts corresponding to zero and positive energy eigenvalues. Although the latter are 
not normahzable, one cannot discard them. Expansion of the eigenfunctions of A over the fixed 
energy states should contain in general an integral over the continuous spectrum. This situation 
differs drastically from the q = 1 algebra case where the Fock space was complete. 

Note that the definition (|5.1(J|) works for > 1 as well since A remains to be the lowering 



operator for positive energy states. Then, the states \a) are normahzable for arbitrary value of a, 
\C{a)r = E,-.{z) = J,. = (-^; q-'U z = ' ' ^ (5.11) 

n=0 \" 1 y Jn ^ 

where Eq-2[z) is another analog of the exponential function, or oV5o('?~^5 —z) basic hypergeometric 
function. In this case coherent states are complete because the Hamiltonian has only a discrete 
spectrum as in . 

For < 1, coherent states formed by the positive energy states should be defined as eigenstates 
of the operator , 

A^a.q) = a\a,q) , (5.12) 

since for E' > 0, not A, but A'^ lowers the energy |T^. Suppose that for some A the states of the 
series ( |5.8| ) are normahzable (this is not so in our case, but in principal it is possible). Then we 
find 

/ oo / n/2 n n{n-l)/2 \\ 

|a, q), = C{a) |0). + ^ -r(- V/^; q')l/'\n)x + ,^ .J -n)A]. (5.13) 



n=l 



A"/2(-z//A;g2)y2 



The normalization constant C{a) is related to the bilateral basic hypergeometric series o^i, 

\C{a)\-^ = oMb; q\ z), h = -z//A, z = -\a\^/\. 
The general bilateral g-hypergeometric series is defined as follows 

Using the Ramanujan sum for the I'j/'i series one can express C[ol) in terms of the infinite 
products 



(g^g')oo(-AgVI« 


2. 2X 


(- 


a 


V^;q')oo 


(-z//A;g2)^(z//|a 


2;g2) 


oo 
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The states (|5.13| ) are normalizable when |ap > u, i.e. a should he outside of the region where the 
states ( |5.10| ) were defined. 

When positive eigenvalues of the abstract Hamiltonian H, H\X) = A|A), form a continuous 
spectrum, the states \a,q) are defined by the integral over |A), 

fo- X\iX,a)\X) Iny^/g^g 

\a,q)= dX d = — - — , (5.14) 

^0 ^(-AgVz/;g2)^ Ing^ 

where xW is an arbitrary function periodic on the logarithmic scale, xIq"^^) = xW- principle 
the continuous spectrum may have an infinity of gaps, each type of gap appearing in the form of 
geometric series. Here we assume that it fills the whole interval < A < oo. From the requirement 
for A and to be hermitian conjugates of each other, with the action (which is defined only up 
to an arbitrary phase factor depending on A) 



A\X) = V^TTX \Xq-^), A^\X) = ^u+Xq^ \Xq^) 
one finds the form of the scalar product 

{X\a) = X6{X-a). (5.15) 



In our case there are two sets of states ( |5.14| ) because the continuous spectrum is doubly degen- 



erate. Expanding the function xW i^ito a Fourier series, we obtain an infinite number of linearly 
independent coherent states of the form 

I v ^ dX X'^' \X) , , 27ris ^ ,^ 

\a,q)s = C{a) 4 = d+-— , s = 0,±l,..., (5.16) 

which are normalizable for |ap > v. The normalization constant is calculated exactly being 



determined by the special case of the Ramanujan g-beta integral |50 



1^/ N|-2 / ciAA^ TT {q '^^]q^)oo ( i^^^^^ 

\C[a)\ - 



Jo {-Xq^ /i>]q^)oo sinvrr (g^;?^)^ \q^ ) 
where r = 2Re d + 1. The fact that there exists an infinity of normalizable states ( |5.16| ) leads to 
an interesting effect in the model of g-oscillator interacting with a classical current [|60|] . 

Consider the structure of coherent states of the first type ( ^.10| ) in the realization ( |5.2| ), (p.3|). 



Denote ipa{x) = {x\a,q). These wave functions are defined by the equation: 

{d/dx + f{x))tlJa{x) = a\f\q\ilja{(ix), (5-17) 
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where f{x) is a solution of ( p.3| ). Unfortunately, a full description of the properties of ipaix) is 
not accessible at present. Nevertheless, several crucial points can be brought into view. 
Suppose < g < 1; then in the a; — > oo limit we get the equation 

ip'c.ix) = ay/qi)a{qx) - ^/utpaix), (5.18) 

which is exactly the pantograph equation encountered in the free particle model. Using the 
corresponding analysis, we conclude that the x ^ oo asymptotics of coherent states is 



-ipaix) = — , K= , (5.19) 

x'^ Inq 

where h{x) is some function satisfying h{qx) = h{x). A similar asymptotics of ipa{x) holds for 
X ^ —oo because in this limit f{x) — > ~\fv. Therefore coherent states are normalizable near 
infinity provided Re k > 1/2, or |a;p < u, consistent with the previous considerations. Due 
to the arbitrariness of h{x,a), asymptotics ( |5.19D corresponds to a countable set of solutions of 
( |5.18|) for fixed a, whereas it is possible to construct only one coherent state by superposing 



discrete spectrum states. Probably it is the requirement for analyticity at x = that determines 
the latter solution uniquely. Then the non-analytical solutions should not be normalizable near 
this point. For Re k < the functions ipa{x) are not bounded, and so not physical. The 
significance of the bounded at infinity but not normalizable solutions of ( p.l7|) , appearing in the 
range < Re k < 1/2, is not yet clear. This problem is related to the analysis of the completeness 
of coherent states which is beyond the scope of the present work. 

Consider the zero modes of the Hamiltonian for the solution of ( |5.3|) determined by the condi- 
tion /(O) = 0. Denoting ipci{x) = {x\cl), we have 

ij'^i{x)=u{x)Mx). (5.20) 

The first several terms of the Taylor expansion of the functions f{x) and u{x) are easily found: 



LUX 



-,2 



— l)uj'^x^ 



fix) = — ^ + — + O(x^), 
^ l + g2 3(l + g2)(i + g4) V 

, , 2u 1u?x^ ^, 4, 

- 1 (1 + g^)^(l + g*) 

Looking for ipci{x) in the form of a Taylor series, we find the first three terms of the odd and even 
wave functions: 

f'cl - 1 + ^4 _ 1 + 3(1 + ^4)(1 _ ^4)2 + 
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« = ^ + W^) ^ 15(1 + - ^ 

Since the operator A maps the space of zero modes onto itself, there should be at least one 

eigenfunction of A. Taking the linear combination of the above solutions we find that there are 
two such functions 

A^ti = ±^v^g-Vcl, (^) = i^:ri^) ± ^v^q-'/'rcf^x). (5.21) 

Evidently are eigenfunctions of the operator as well, A^^p"^ = ^iq^/vilj^i. The eigenvalues 
of A and A'^ are not complex conjugate to each other. Hence, the "condensate" representations of 
the g-oscillator algebra do not belong to the discrete spectrum, actually, ^/'^ are not even bounded. 
This can be verified by the direct estimation of the asymptotics of these functions using the fact 



that they satisfy equation (|5.17| ) with a = ±2y^/g. We, however, choose a slightly different 



approach. Consider the eigenvalue equation for A"^ 

{-d/dx + f{x))^^fi{qx) = TtqV^^fiix). (5.22) 
Eliminating the derivative part from ( |5.17| ) and ( [5. 22] ), we get: 



(fix) + q-'f{q-'x))ijti{x) = ±Zv^g-i/2 (^,^(gx) - ^/jfiiq-'x)) , (5.23) 

i.e. zero energy wave functions satisfy the second order purely finite-difference equation. Since 
for X — > oo, /(x) ^ -0/, the leading asymptotics of ip"^ are found as solutions of the free finite- 
difference Schrodinger equation, which are 

ij^ix) = x^(x)e^'=±i'^^ x^{qx) = X^{x). (5.24) 

Substituting this Ansatz into (|5.18|) with a = iiy^/g, one finds 2ik± = 1 =i= ivr/ Ing. The derived 
asymptotics of ipfi{x) look like Bloch wave functions for a particle in a periodic potential with 
coordinate variable being Inx. Note that the quasimomenta k± are complex, which forces the 
wave function to increase at infinity as ^/x, unlike the u{x) = case when a symmetric constant 
wave function was bounded. A more detailed consideration of such non-standard implementation 
of the Bloch's theorem is given in section 8. 

Let us discuss briefly the — 1 < g < case. Due to the parity symmetry one has the relation: 
^^lg<o ~ ~^'^\q>o- means that the g > coherent states \ia, q) and |— ia, q) provide independent 
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eigenstates of A^I^^q. The parity transformation changes only the sign of a. Therefore, the 
eigenstates of A for g < are given by the parity coherent states ( |2.6| ) where on the r.h.s. one 
has g-coherent states for g > 0. Moreover, the precise reahzation ( |2.2| ) and superposition ( p^.6D 
appear in the g — ^ — 1 hmit of the g-oscillator system under consideration [13]. Indeed, the general 



solution of equation ( |5.3| ) for g = — 1 is /(x) = u;x/2, and the variation of g from 1 to —1 performs 
a transition from the canonical coherent states to the parity coherent states (|2.6|) . Note however, 
that if one takes solution of ( ^.3|) satisfying asymmetric initial condition /(O) 7^ 0, then the g ^ — 1 
limit simply does not exist and for — 1 < g < the parity operator does not help in the analysis. 
Coherent states of the second type are defined as eigenstates of the operator ( ^.12[) : 

{—djdx + /(a;))yg-?/'Q,(gx) = ail)a{x)^ (5.25) 



where we use the same notations as in ( |5.17| ). In the x ^ 00 limit one gets again the pantograph 



equation, but with the scaling parameter g ^: 

i)'^{x) = -aq-^/^iJa{q'^x) + q-^^i){x). (5.26) 
From the preceding section we know that this equation has solutions with asymptotics 



^ ^ M:^^^ ^^^^^ ^ ^ ^ In v^/a ^ 

x'^ In g 

These coherent states are normalizable provided Re k > 1/2, or \a\'^ > v. Their expansion over 
the continuous spectrum states was presented in the abstract form ( p.l(j| ). For Re k < 0, or 
lap < qv the functions ipai^) are not bounded. For < Re k < 1/2 we get again wave functions 
that are bounded at |a;| — ^ 00 but unnormalizable. It looks like they do not have an expansion over 
the eigenstates of the Hamiltonian and the reason for this needs clarification. In the considered 
models there are no normalizable coherent states for the circle |ap = z/. In principle it is possible 
that zero modes of a Hamiltonian H are normalizable, in which case it is natural to count them 
among the coherent states (otherwise there will be no completeness). 
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6. General class of self-similar potentials and their coher- 
ent states 



It is well known that the one-dimensional Schrodinger equation 

Hip{x) = —ip"{x) + u{x)ip{x) = Xip{x), (6.1) 

has an important non-quantum mechanical application in the theory of non-linear evolution equa- 
tions. In particular, the Korteweg-de Vries (KdV) equation can be solved with the help of the 
inverse scattering method for two classes of initial conditions u{x,t = 0). The first one consists 
of the potentials u{x) satisfying the restriction + |a;|)|u(a;)|(ia; < oo, which guarantees that 

the number of bound states is finite. Reflectionless potentials with N discrete eigenvalues are the 
simplest examples from this family. Since they generate A^-soliton solutions of the KdV equation 
they are called the soliton potentials. The second class is related to non-singular periodic (or 
quasiperiodic) functions, u{x + I) = u{x), characterized by the presence of gaps of finite width 
in the spectrum of (|6.1|) . These finite-gap (hyperelliptic) potentials are reduced to the solitonic 



ones in the limit I —>■ oo. They can be thought of as superpositions of an infinite number of solitons 
( "periodic solitons"), but there is no scattering problem for such objects, i.e. the solitary character 
of the ingredient waves is lost. There is a third relatively simple class of self-similar solutions of 
the KdV equation related to the Painleve transcendents requiring a different treatment. 

Recently, reflectionless potentials with an infinite number of discrete levels have been system- 
atically considered in |6^, ||, |, |13|, |5^. For x — > oo such potentials decrease slowly and 



the standard inverse scattering method does not produce constructive results. These potentials 
deserve to be named the infinite soliton ones, since they do not reflect and may be approximated 
with some accuracy by the A^-soliton potentials (A^ < oo). Moreover this class absorbs the finite- 
gap potentials which emerge for special limiting values of parameters. The related problem of the 
approximation of confining and band spectrum potentials with the help of reflectionless potentials 
was discussed earlier in |67| . 

Let us describe briefly the factorization method 0, ^] that allows us to find the particular 
subclass of infinite-soliton potentials characterized by the g-deformed symmetry algebras. This 
method was invented in quantum mechanics by Schrodinger, it is deeply related to the Darboux 
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(Backlund, dressing, etc) transformations for linear differential equations. Within this approach, 
one takes a set of Hamiltonians, 

Lj = -d^/dx"^ + Uj{x), j = 0,±1,±2,... (6.2) 

and represents them as products of the first-order differential operators. 



/i; = -^ + /,(x), Aj = ^ + m^ (6-3) 



up to some constants A^: 

L, = AjA- + \j, (6.4) 
i.e. Uj{x) = fj{x) — fj{x) + \j. Then one imposes the following intertwining relations: 

= A+Lj+i, AjLj = Lj+iAj, (6.5) 

which constrain the difference in spectral properties of Lj and Lj+i and are equivalent to the 
equations: 

Aj_,,Aj^, + A,+i = AjAj + A,. (6.6) 
The same results are obtained if one starts from the linear equations 

Lj-^j = ipj+i = Aj-^j. (6.7) 

The compatibility condition of ( |6.7|) is given by ( |6.5|) . Resolving the latter relations, one finds 
Xj as the integration constants such that the factorization ( |6.4| ) takes place. These are the basic 
ingredients of the factorization method which allows us to construct new solvable Schrodinger 
equations from a given one. 

Substitution of (|6.3|) into (|6.6|) yields the chain of differential equations: 



fj{x) + fj+i{x) + ff{x) - ff+^{x) = fij, fij = Xj+i~Xj, (6.8) 

which is called the dressing chain |^, Any spectral problem with a known nontrivial discrete 
spectrum generates some solution of ( |6l8| ) such that the ordered discrete eigenvalues are given 
by the constants Xj. The factorization method works with the inverse problem - to find such 



solutions of ( |6.6|) or ( |6.8|) for which the spectrum of the associated Schrodinger operator will be 
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determined automatically. Note that in some cases the spectrum can be found even if Xj do not 
belong to it, e.g., such a situation can take place for hyperelliptic potentials [^, p9|-|[7l[]. 



The variable j was playing above the role of a label. It could therefore be removed in favor 
of other notations, e.g. fj = f, fj+i = f, etc. However, it is convenient to think of j as a 
discrete set of points on a continuous manifold. Then one can consider j (or its function) as a 
continuous variable and look for solutions of (|6.8| ) in series form fj{x) = J2 9k{x)j'' ■ Infeld and 
Hull [§] have considered such an Ansatz and have found that the series contains a finite number of 
terms iff fj{x) = a{x)j + l3{x) + 7(a;)/j, where a,/3,7 are some elementary functions. This does 
not mean that the infinite-series solutions are meaningless, it indicates rather that truncation of 
the series is related to some simple symmetries 0. It is difficult to work with formal power series, 
often even their convergence is not known. Analysis of solutions characterized by separation of 
variables or by some special dependence on them is essentially simplified. In the Lie theory of 
differential equations solutions of such type are called the similarity, or self-similarity solutions. 
Unfortunately for the differential-difference equations there is no complete theory of such solutions; 



a subclass of them can be found using the methods developed for purely differential equations []72 
An additional complication associated with equation (|6.8|) consists of the fact that there are two 
unknowns, fj{x) and fij, i.e. the system is highly underdetermined. 

According to definition, (self-)similarity solutions are the solutions invariant under symmetry 
transformations of a given equation. The potentials we are interested in appear as fixed points of 
the combination of an affine transformation of the coordinate, x —>■ qx + I, and a shift along the 
discrete lattice, j —>■ j + N. Indeed, the change in numeration of solutions by an integer, maps 
solutions of ( |6.8| ) to the solutions, fj{x) — > fj^j\f{x), fij — > fij+N- The same is true for the affine 
group, fj{x) — > qfj{qx + I), fij — ^ q'^fJ'j- We may look for the class of solutions invariant under 
both these symmetry transformations: 

fj+Nix) = qfj{qx + l), Hj+N = q'^IJ'j- (6.9) 

These relations define a class of potentials that we are going to analyze below. The simplest 
reduction of such type, fj{x) = q^ f{q^x), Xj = q"^^ , corresponding to = 1, / = in ( p.9| ), 
has been found by Shabat p. Actually, the general class of closures of the dressing chain ( |6.9| ) 
has been introduced in a way different from the above [0, namely, from g-deformation of the 



parasupersymmetric quantum mechanics based upon some polynomial algebras [^ . 
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At the operator level, the relations ( |6.9| ) lead to the Schrodinger operators with non-trivial 
g-deformed symmetry algebras. Let us consider the products: 

m; = . . . M- = AT+^.i . . . A-^,A-, (6.10) 

which generate the intertwinings 

LjM+ = M^Lj+N, MjLj = Lj+nM-. (6.11) 



The structure relations complimentary to ( |6.11|) appear as 



Af--1 N-1 



M+M7 = n (L, - A,+fc), M-M; = n (W - (6-12) 

k=0 k=0 

Equations ( |6.11| ), ( |6.12| ) can be rewritten as the higher order polynomial supersymmetry algebra 
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N-l 2 

{Q^, Q-}= liiK-Xk), (g^) = [K, Q^] = 0, (6.13) 

k=0 

where 

\0 ) VMq- 0/ V Ln, 

the index j = was fixed for simplicity of notations. Presence of nonlinearity brings essentially 
new features with respect to the standard supersymmetric quantum mechanics case corresponding 
to = 1. 

The identities (|6.11|) show that if the operators Lj and Lj^^ are related to each other through 
some simple similarity transformation, e.g. 

Lj+N = q^TL^T-' +00, (6.14) 

where T is some invertable operator, then the combinations 

B+ = M+T, B-=T-'M-, 

map eigenf unctions of Lj onto themselves, i.e. they describe symmetries of Lj. The form of T 
is restricted by the requirement for the L^'s to be of the Schrodinger form. The closure ( |6.9|) 
corresponds to the choice of T as the affine transformation operator, T f[x) = J\q\f{qx + I). 
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Fixing the indices and removing their irrelevant part {L = Lq, = -B^), we get the symmetry 
algebra H, 0] : 

LB+ -q^B+L = ujB^, B' L - q^LB~ = tuB' , (6.15) 

N-l N-l 

B+B- = n - >^k), B~B+ = n {q^L + UJ- (6.16) 

fc=0 fc=0 

After the shift of the zero energy point, 

H = L — Ek = Afc — 



l-g2' - i-q2 

the algebra takes a simpler form 

HB^ = q^^B^H, 

N-l 

B+B- = HiH-Ek), (6.17) 

k=0 
N-l 

B-B+ = Y[iq^H-Ek). 

k=0 

This algebra was already met in the discussion of coherent states of the free nonrelativistic particle. 

Let us write out explicitly the system of nonlinear differential equations with deviating argu- 
ment that one needs to solve in order to find the explicit form of the self-similar potentials: 

^ ifoix) + fi{x)) + /o (x) - f^{x) = /io, 

^(/1(X) + /2(X)) + A^(X)-/|(X) = /ii, 



^{fN-i{x) + qfo{qx + l)) + fN_,{x)-q'f^{qx + l) = (6.18) 

Note that the limit g ^ 1 is not trivial. For non-zero parameter I we get a realization of the 
algebra (|6.15|) , ( |6.16| ) at g = 1 which generalizes the one described in |^8[. Below we assume 



that / = (for q ^ I this corresponds to the fixed point reference frame but for q = I this 
gives a nontrivial simplification). The zero potential model considered in section 4 corresponds 
to a particular solution of this system, fi{x) = [3j = const, which determines thus the simplest 
self-similar potential. 

Consider some examples. The = 1 case describes a g-deformation of the harmonic oscillator 
potential, since for g = 1 one has f{x) = ^x/2 and u{x) oc x^. The N = 2, g = 1 system coincides 
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with the conformal quantum mechanical model [ffSf, 



fo,i{x) = ±- + (3x, 7 = 777 \ V' ^= — ~A ' ^-^9 

X 2{no + Hi) 4 

Uo,iix) = + p^x^ _ ^(1 ^ 27) + Ao,i. 

This is the singular oscillator potential whose physical spectrum consists formally of two arithmetic 
series but only one of them is physical due to the boundary condition at x = 0. Coherent states 
of this model defined as eigenstates of the symmetry operator B~ were constructed by Barut and 
Girardello (for some amendments, see |7^). The N = 2, q = —1 system obeying the same 



symmetry algebra will be considered in the next section. 

It is natural to call the physics of the N = 2, ^ \ models as the g-deformed conformal 
quantum mechanics because their symmetry algebra is 1). A general solution of the basic 

equations (|6.18| ) is not available in a closed form already for = 1. Let us find /o,i(a;) as formal 
series near x = 0. Consider first singular solutions with the pole type singularity. Such solutions 
appear to be odd functions with a simple pole at zero: 

/o(x) = - + f:&.x^^-\ = - - + f:c.x^^-\ 

i=l i=\ 

Pi 2z - 1 + 2a ' ^ ''^ + - A. 2z - 1 - 2a ' ^ ^ ^ 

where z = 2, 3, . . . , and a is an arbitrary parameter; the coefficients 6i, Ci have the form 



l-g2Vl + 2a l-2ay l-g2\^l-2a 1 + 2a / 

In general, the series diverge at g — 1, but for the special choice of a one gets the truncated 
solution (|6.19| ). In the limit q 0, the function qfo{qx) does not vanish, qfo{qx) — > a/x, and then 
the potentials uq.i are obtained by simple dressing of the potential a(a + l)/x^. It is natural to 
expect that for < g < 1 there exist yUo,i such that the series converge for arbitrarily large x. The 
condition a(a + 1) > 3/4 guarantees that normalizable wave functions and their first derivatives 



vanish at zero The spectrum of such a system would arise from only one geometric series 
(for the same reason that there is only one arithmetic series for the singular oscillator). 

Due to the presence of the singularity one can restrict the space to be a half-line, < x < oo, 
and interpret x as a radial coordinate appearing from the separation of variables in a three- 
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dimensional Schrodinger equation. The constant a aquires then an interpretation of orbital mo- 
mentum of a particle when it takes integer values a = I. Then in addition to the standard term 
/(/ + there is a very complicated dependence of the potential on the quantum number /, 

indicating that in this picture the "g-deformed Laplacian" is a nonlocal operator. 

For the solutions that are non-singular at zero one has /o.i = Z^i^o ^i"'^"*^*? where &o°'^'* are 
two arbitrary constants. Again, in general the series diverge at g ^ 1. A particular choice of 
initial conditions gives the solution which in this limit corresponds to (|6.19|) with a coordinate 



shift. In the < g < 1 region, the solution that is nonsingular for all real x defines an infinite 

soliton potential whose spectrum is composed from two independent lowest weight irreducible 

representations of s-Ug(l, 1). In the limit g ^ it shrinks to the smooth two-soliton potential. It 

is this solution that reduces to the g-oscillator one ( |6.12| )-( |6rT6D (with g replaced by g^^^) after 

the restrictions fi{x) = q^^"^ fo{q^^'^x), jii = qjiQ. Probably without such a restriction this solution 

does not have a g — ^ 1 limit. In any case, at g 1 the spectral series become equidistant, which 

means that the potentials start to be unbounded at infinity. Obviously a more complete and 

rigorous analysis of the structure of the solutions of the N = 2 equations is necessary. 

When g = 1, already the = 3 case corresponds to transcendental potentials, namely, the 

fj{x) depend now on solutions of the Painleve-IV (PIV) equation [ |6^ : 

1 

fo{x) = -ux + f{x), u = fiQ + fii + fi2, (6.21) 

f" = ^ + If + 2^^/' + (l^'^' + f^2-f^o)f- §j. (6.22) 
Thus, a simple generalization of the harmonic oscillator, characterized by the split of its linear 
spectrum onto A^ independent terms, is connected with highly nontrivial ordinary differential 
equations whose solutions are transcendental over the solutions of linear differential equations 
with coefficients given by the rational and algebraic functions. For the A^ = 3, g 7^ 1 system it is 
not possible to reduce the order of the equation (i.e. no first integral is known). Because in this 
case there are solutions of ( |6.18|) reducing in the limit g — 1 to the PIV functions, it is natural to 



refer to this system as to the g-deformed PIV equation. 

The notion of g-periodicity (|6.18|) and the corresponding algebraic relations are central in this 
section. Suppose that superpotentials fj{x) are not singular and that the operator is well 
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defined (-B^ = [B^Y) and has normalizable zero modes (the necessary condition for this is 

B-\l)=0, {l\m) = Sim, l,m = 0,l,...,N -1. (6.23) 

Then |/) represent the first bound states of self-similar potentials and the whole discrete spec- 
trum consists of N independent geometric series: 

H\n) = En\n), En+N = q'^En, En < En+1 < 0, 

EkN+i = Eiq'"', / = 0,l,...,iV-l, fc = 0,l,...,cx), (6.24) 
\Nk + l) = Ck{B+)''\l), {Nk + l\Nk' + /') = 5kk'5iv. 
The normalization constants Ck are determined up to an arbitrary energy-dependent phase factor 

N-l k 

= i^l[ n (1 - q'^'Ei/E,), V = {-ifEoE, . . . > 0, 

s=0 m=l 

and similar freedom exists in the action of raising and lowering operators 

N-l 



B+\kN + = n \/Eiq^(-^+^^ - E,\{k + l)N + /), 

s=0 



N-l 

B-\kN + l)= W ^Eiq^^ - E,\{k-l)N + 1). (6.25) 

s=0 

In the "crystal" limit g — > there remain only first levels and we get the general A^-soliton 
potentials. 

There may be intermediate situations, when only some of the zero modes of B~ satisfy neces- 
sary boundary conditions (the singular oscillator model is a known example). Then some of the 
spectral terms disappear; we do not consider such possibilities in detail here. Note that for g > 1 
the above formulas probably only have formal meaning in the differential Schrodinger equation 
case due to a too rapid growth of spectrum. Analysis of the A^ = 1 potential showed that, for g > 1, 
it has singularities, and thus the symmetry operators are not well defined Q. The finite-difference 
Schrodinger operators, or the Jacobi matrices have a more rich spectral structure |37, 78, 59 1. In 



particular, the spectrum may consist of a discrete set of points only and be compact; it may grow 
exponentially fast, and the geometric series may accumulate near zero both from below and from 
above, etc. 
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Let us summarize what one can do with the help of Darboux transformations. Starting from 
the zero potential u{x) = 0, and performing dressing transformations, one can get A^-soliton 
potentials. Starting from the latter and performing an infinite number of dressing transformations 
in the particular self-similar manner, we get a subclass of infinite soliton systems. Taking the 
g — > 1 limit we derive the Painleve-type potentials, some of which degenerate into the finite-gap 
potentials in the limit uj ^ 0. All these complicated systems are thus equivalent in some sense to 
the free Schrodinger equation. 

Consider coherent states of the algebra (|6.17| ) defined as eigenstates of the operator , 



B-\a) =a'^\a), (6.26) 

where for convenience we denote the eigenvalue as the A^th power of a. Since B~ has A^ linearly 
independent zero modes, there are A^ independent states: 

fc=0 11^=0 llm=ll-^«? - ^s) ' 

|^|2(fciV+0 



fc=0 lls=0 llm=ll-£^«'? i^s) 

\ \2l / 0, . . . , 2 1 I |2A^ / 

= |a| ArV?Ar-l ; ^ ; ? , 2: , Z = \ol\ jv. 
• • • , Oat-i / 

h\ = Eiq^Eo, . . . , = EiqyEi^,, . . . , b'^_, = Eiq^E^^i. 

These coherent states are normalizable for < z/. In order to find coherent states corresponding 
to the limit g 1 it is necessary to replace the g-products by the Pochhammer symbols. The 
spectrum of the resulting system is found from the formula En+N = En + uj. The hierarchy of 
potentials with such a linear spectrum was investigated in |^ . Representation theory of the g = 1 



algebras was considered in |7S|, Their possible physical applications were discussed also in 



8T| , p^ , p3[ . In the limit Ei Eq + Icu/N, I = 0,1, . . . , N — 1, one finds the harmonic oscillator 
such that the states \ai) constructed above coincide with the generalizations of even and odd 
coherent states (|2.14|) . 

The positive energy states are described by the following representation 

H\X) = A |A), A>0, 

N-l 



B^\X) = n VV-i5;.|Ag2), 
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N-1 



|A) = n ^^-Es\M-^). (6.28) 

s=0 

When the states generated by from |A) for some fixed A are the only eigenstates of the 
Hamiltonian, generahzation of the coherent states (|5.13|) has the following form: 



«)a = C{a) E J, „ |Ag-^") (6.29) 



^nN qNn{n-l)/2 

oo ra/2 N-1 ^ 

n=l " s=0 / 



where 



oo / 'i\Nn„Nn(n—l)^n oo n N—1 

n=0 lls=0 \^sl^,q )n n=l l"l s=0 

= oV^7v(&o, . . . , bN-i; q\ z), z = (-l)^|a|2^/A^, h = Ei/X. 

These states are normalizable for > z/. 

Analogously, one can find the countable set of coherent states in the case when positive energies 
form the continuous spectrum: 

where the normalization constant C{a) is given in ( 4.20 ) with jSf = —Ei. Zero modes of the 
Hamiltonian H may be eigenstates of the operators A and in the same way as in the g-Weyl 
algebra case; there is no need to describe them again. 

We conclude that coherent states of the self-similar potentials are related to the basic hyper- 
geometric series nVn-i{,z), qiPn{z) or to some integrals over the latter functions. However, we do 
not know the explicit form of these coherent states because the Hamiltonian eigenstates involve 
new complicated transcendental functions whose complete analytical structure is not accessible 
at present. Let us mention in passing that the finite-difference analogs of some of the above 
self-similar potentials associated with the discrete Schrodinger equation have been described in 
| 78|| . The discrete dressing chain in this case happens to define a discrete time Toda lattice [|59[ . 



A g-deformed supersymmetric interpretation of the models considered in this and the preceeding 
two sections was given in p4|, 0, |74| . 
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7. Two particular examples 



Consider in more detail the N = 2, q = —1, and a special subcase of the = 3, q = 1 closures 
of the dressing chain. In the first case, the system of equations determining superpotentials has 
the form: 

^(/o(x) + + f^{x) - fl{x) = /xo, (7.1) 

- foi-x)) + fl{x) - fl{-x) = /.I, 

where we assume that /iq + /Ui 7^ 0. From the algebraic point of view this case is related to the 
su{l, 1) algebra which serves as the formal spectrum generating algebra. If one assumes that /o(x) 
(or /i(a;)) is antisymmetric, /o(— a;) = —fo{x), then our system is equivalent to the N = 2, q = I 
case, or to the singular harmonic oscillator model (|6.19|) . The coordinate region of this problem 
is restricted to the half-axis which forbids the parity operator. Indeed, the wave functions are 
not single-valued near x = for noninteger values of the parameter 7 and the action of the 
parity operator is not well defined. As a result, the eigenvalue equation for the lowering operator 

= P{d? /dx'^ + . . .) is also not well defined (for integer 7, half of the wave functions are singular, 
the coherent states have fixed parity, and the action of and B^ P may differ only by sign). 

It is easy to see that fo{x) or fi{x) cannot be symmetric, therefore the non-trivial solutions of 
( [7. 1|) do not have fixed parity. In spite of its simplicity, the system ( [7. 1|) is hard to solve. Let us 
represent fj{x) as sums of symmetric and antisymmetric parts: 

/s,a(x) = ^(/o(a;) ± /o(-x)), ^?,,,(x) = ^(/i(x) ± h{-x)\ (7.2) 



and substitute this splitting into (7J). Then it is not difficult to find from these equations and 
their x —>■ —x partners two "integrals" 

fa{x) + ga{x) = ax, o- = ^(/io + /ii), (7.3) 



and 



fsix) + faix) = gl{x) + gl{x) + r, r = l(/io - /ii). (7.4) 



As a result, it is possible to express fa and Qa via their symmetric partners 

fa{x) = ^^{gl-f! + a'x' + r 
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9a{x) = ^{f^-9'.+a'x^-r). (7.5) 
Eventually the original system of functional-differential equations is reduced to the form: 

fiix) = 2gs{x)g,{x), g^x) = -2f,{x)Ux). (7.6) 

This system of ordinary differential equations is related to a Painleve-V equation, which follows 
from the fact that we are considering a subcase of the = 4, q = 1 closure analyzed in [Q. 

For the special choice of initial conditions /a(0) = 5'a(0) = one has (0) = g'^{0)+T. Then the 
singularity at x = cancels and one obtains a formal nonsingular solution of the equations ( [7.6|) . 
It is possible that this solution defines a smooth potential growing indefinitely when —>■ oo. 
Under this asumption, consequences of the presence of the symmetry algebra sm(1, 1) in this model 
are different from those for the singular oscillator potential. Namely, the spectrum consists now of 
two arithmetic series, each being determined by irreducible representations of the su{l,l) algebra. 
The coherent states \ai), / = 0, 1, or even and odd coherent states, are both physical but they are 
not eigenstates of the parity operator because the potential is not symmetric. 

Another example that we would like to discuss in more detail concerns the = 3, g = 1 closure 
considered in [^] (note that in this paper all Hamiltonians L, admitting A^th order differential 
operator 5 as a symmetry operator satisfying the relation [L,B] = B, were characterized). The 
PIV function appearing in this context has an infinite set of rational solutions emerging for the 
specific choices of the parameters /xq, /^i, |@- Consider the following rational solution: 

fiix) = (7.7) 
2 X x"^ + 1 

which corresponds to the constants /xq = 4, /xi = 1, /i2 = —2. The symmetry algebra has thus 
the form 

[L,5±] = ±35^, 

fi+fi- = L(L-4)(L-5), (7.8) 
B-B+ = (L + 3)(L- l)(L-2), 
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where the Hamiltonian L is 

(PA 8 1 2 3 , 

+ 1 (a;2 + 1)2 4 2 

There are three normahzable eigenstates of the lowering operator B~ corresponding to the lowest 
energy state, and second and third excited levels. Acting by upon these three "vacua" one 
finds the whole spectrum. It consists of three arithmetic series forming the sequence 0, 3, 4, 5,6,.. ., 
i.e. there is a hole between and 3 after which the spectrum becomes equidistant (the first excited 
state with the energy 3 is generated by B^: |1) oc -B+|0)). 

The explicit form of the coherent states, B^\a) = a^\a), is determined from the equation 



2 



X"'i^) - ^^X"{x) + ^^2^^y X{x) = a'xix), (7.10) 

where 

X{x) = (x2 + l)e"'/^(x|«). 

The author does not know any simple solution of this equation. Each of the three linearly inde- 
pendent solutions defines a particular coherent state. Evidently, they have some common origin 
with the cubic root-of-unity superpositions of the coherent states (p.l4|) . 



The structure of spectrum hints that there are some additional symmetries in this model, and, 
indeed, there is a different set of raising and lowering operators satisfying the different algebra 
found earlier in [Q. It corresponds to a different solution of the dressing chain for the same 
Hamiltonian ( |7.9| ): 



fo{x) = -Mx) = - + ^:^, fi{x) = ^, (7.11) 



for which /zq = 3, /xi = — 2, fj,2 = 0. Now the symmetry algebra has the form 



[L,5±]=±S±, 

B+B- = L{L-3){L-1), (7.12) 

B-B+ = (L + 1)(L-2)L. 

The operator B~ has only two normahzable zero modes corresponding to the first two levels of the 
Hamiltonian. The principally important feature of this realization is that the lowest energy state 
is also a zero mode of the raising operator. As a result, it is the first excited level which serves 
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as the lowest weight vector of the discrete series representation of the above algebra from which 
the raising operator creates an infinite tower of states. In this picture the spectrum consists 
of only one arithmetic series with step 1, and the isolated zero energy state forms a trivial one- 
dimensional irreducible representation of the algebra. One may thus conclude that the spectrum 
generating algebra of a given potential, and hence the coherent states, may be non-unique. 

The explicit form of the coherent states of the second symmetry algebra are determined by the 
equation 

X'"ix) - a: (l + ^) x"{x) + (2 + ^ - J^^) x{x) = (7.13) 
where x{x) is defined in the same way as in the first case. Simple solutions of this equation are also 



not known to the author [^. There should be only one physically acceptable solution because 
now the coherent states have a unique expansion over the Hamiltonian eigenstates for a ^ (i.e. 
there is no connection with the root-of-unity superpositions). An analogous situation has been 



described recently for a different particular Hamiltonian with equidistant spectrum |88]. It would 



be interesting to find physical characteristics of the coherent states ( |7.10| ), ( |7.13| ) and similar ones 



8. Schrodinger operators with discrete scaling symmetry 



Spectral theory of the one-dimensional Schrodinger equation ( |6.1|) for the bounded periodic po- 
tentials, 

u{x + l)=u{x), (8.1) 



is well known |^. The Bloch, or Floquet theorem states that there is at least one wave function 
tplx) which can be represented in the form: 

^(x) = e^'^Xx), (8.2) 

where (f){x) is a periodic function, 0(x + I) = (pix). When iIj{x) is a bounded function it defines a 
physically acceptable state. This requirement restricts the quasimomentum k to be real. 
Consider a class of Schrodinger operators, which satisfy the relations 

T^H = q^HT\ HT = q^TH, T^T = TT^ = 1, (8.3) 
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where T is the unitary fixed affine transformation operator ( [4 .91) . The algebra ( p.3| ) corresponds to 
the degenerate case of ( |6.17] ) when the polynomials of H on the r.h.s. are replaced by a constant. 



It is equivalent to the following constraint upon the potential: 

q^u{qx + 1) = u{x), (8.4) 

which for g = 1 is the periodicity condition ( |8.1| ). Note however that the limit g ^ 1 is not defined 
uniquely. After the shift u{x) — > u{x) + cr/(l — g^) the condition ( ^.4|) is reduced in this limit to 
u{x + 1) = u{x) + a, which shows that u{x) is the sum of a periodic and Airy potentials. Only the 
additional requirement a — > leads to the standard Bloch-Floquet theory. 

For the shifted potential v{x) = u{x + x*), where x* is the fixed point, qx* + I = x*, the 
constraint ( p.4| ) is converted into q^v{qx) = v{x). This means that 

v{x) = h{x)/x'^, h{qx) = h{x), (8.5) 

where h{x) is an arbitrary function obeying the indicated property of periodicity on the logarithmic 
scale. In the following we restrict ourselves to the case when h{x) is bounded for < x < oo, 
which is a natural coordinate region of the problem. In this case the potential v{x) is singular 
at a; = and vanishes in the limit x —>■ oo. The value of h{x) is not defined at zero unless it 
is a constant. This means that one cannot soften the singularity by requiring h{x) to vanish at 
a; = 0. The formula ( |8.5| ) defines a class of potentials whose spectral theory we propose to call the 
g-Floquet theory. 

It is instructive to rewrite the resulting Schrodinger equation using the change of variables 
X = expy and renormalization of the wave function iIj{x) = \^x{y)- 

- + {h{y) + \)x{y) = Ae^^x(2/), (8.6) 

where h{y + Ing) = h{y) is now a periodic function. This equation has a form similar to the 
original one but the normalization condition for bound state wave functions xiy)^ 







ix)\'dx= e'y\x{y)\'dy = l, (8.7) 

J —oo 

contains a nontrivial factor e^^ in the measure, which diverges when y oo. 
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The presence of the non-trivial symmetry constrains the structure of the wave functions. The 
operator T acts upon the wave functions as follows 

2 

Tipi{x, A) = y/qi)i{qx + /, A) = ^ Cij^jj^x, g^A), (8.8) 

i=i 

where ilJi{x,X) are two linearly independent solutions of the Schrodinger equation. On the one 
hand, this g-difference equation shows that if ipi^^ ^) is normalizable for some fixed A = Aq, then 
T does not change its normalizability. On the other hand, the resulting wave function has the 
eigenvalue A = g^Ao, which thus also belongs to the physical spectrum. Applying the operator 
one gets a wave function with the eigenvalue A = q^^Xo- This simple argument shows that if our 
Hamiltonian has physical states with negative energy, then it exhibits the "fall onto the center" 
phenomenon [^, the energy is not bounded from below. 

The class (|8.5| ) unifies periodic potentials with the conformal one corresponding to the special 
choice h{x) = const Indeed, the Hamiltonian H = —(P/dx"^ + h/x^, enters the following 

conformal symmetry agebra: 

[H, D] = iH, [K, D] = -iK, [H, K] = 2tD, (8.9) 

where D = i{x,d/dx}/4 and K = are the hermitian scaling and special conformal transfor- 
mation generators. In this case dilatation by arbitrary parameter maps solutions of the Schrodinger 
equation onto the solutions, which means formally that the spectrum of H is purely continuous 
(such a situation holds only for h > —1/A). If h{x) ^ const the discrete spectrum consists of a 
number of geometric series infinite in both directions, with A = and A = — cxd as the accumu- 
lation points. Note that ( ^.51) , with some unknown h{x), was suggested as the asymptotic form 



of the self-similar potential with one geometric series of bound states truncated from below 

In the more general case one has N such series it is a matter of conjecture that in all cases 

asymptotics of the potentials have a similar form. 

From the symmetry point of view we have the following situation. The simple conformal 
potential has a large group of continuous symmetries which fixes the dynamics completely. We 
can add an interaction term which does not remove all symmetries, but preserves a discrete part 
of them. This situation is similar to the case where a free particle is put into a periodic potential 
- instead of the group of continuous translations one has only a discrete subgroup of it. Note that 
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the generator of the discrete subgroup is now considered as an element of the symmetry algebra 
characterizing the spectrum - this is one of the ways of building quantum algebras out of Lie 
algebras. 

Due to the non-analyticity of the potential at x = 0, the accurate spectral analysis requires 
some rigorous mathematical tools. The main problem is to find restrictions upon h{x) for which 
the Hamiltonian is self-adjoint, or for which it has self-adjoint extensions. Here we give only a 
qualitative picture. From ( ^.5|) we see that there are three essentially different regions of the 
spectrum. Positive energy states form continuous spectrum because the potential goes to zero 
sufficiently fast. Degeneracy of these states depends on the boundary condition imposed upon the 
wave functions at x = 0. Only for boundary conditions that are invariant under the taken scaling 
transformation, does the operator T represent a physical symmetry. 

The second region of energy corresponds to zero modes of the Hamiltonian. It is especially 
interesting because upon these zero modes T commutes formally with the Hamiltonian. For A = 0, 
( ^.6[ ) looks like the standard Schrodinger equation for a periodic potential with spectral parameter 
equal to —1/4: 

-X"{y) + h{y)x{y) = - Kv + in g) = h{y). 

If both independent solutions of this equations are bounded, i.e. if the eigenvalue —1/4 belongs 
to the permitted band of h[y)i then the original wave function does not describe physical 
states, being unbounded at infinity. For boundedness of il}[x) at x = oo the quasimomentum k 
of the Bloch eigenfunction xW) = + Ing) = has to have an imaginary part Im 

k > 1/2 (for Im k = 1/2 one has iplqx) = q'^^'^^il){x)). Normalizability of il){x) near x = imposes 
the restriction Im A; < 1. For Im A; = 1, wave functions il){x) are generalized eigenfunctions of T 
with eigenvalues of modulus 1, as it should be for a unitary operator. The zero energy states of 
the Hamiltonian H may thus belong to the continuous spectrum when the quasimomentum k of 
the Bloch function x{y) is complex and 1/2 < Im /c < 1. 

For any A < there is a wave function il){x) which is normalizable at infinity, so that it is 
the boundary condition at x = which determines the quantization of the spectral parameter. 
If h{x) > then there are no bound states at all. When squares of the absolute values of both 
solutions of the Schrodinger equation are integrable near zero, the Hamiltonian is not self-adjoint. 
Its self-adjoint extensions are fixed by the requirement of a particular dependence of wave functions 
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near zero. 

Let us consider the case h{x) = const, which exhibits already basic features of the g-Floquet 
theory. This model was considered by Case |^ who has shown that when h < —1/4 one can choose 



an orthonormal basis of states such that the energies of bound states form an infinite geometric 
series in both directions. Such a spectrum is stipulated by the presence in the model of a discrete 
scaling symmetry. More precisely, the operator D in (|8.9|) is not a physical symmetry operator any 
more, but the group element g^*^ is - it maps wave functions of the discrete spectrum onto each 
other. Let us consider in more detail how this situation appears using the Wronskian technique 
of self-adjoint extensions of singular Schrodinger operators |9^. According to this approach. 



one takes two solutions of the original equation for some fixed eigenvalue satisfying 0i = cos a, 
(f)[ = sin a and 02 = sin a, 02 = — cos a at some regular point x = 6 (evidently, 0i02 — 0'i02 = — 1). 
Then it is necessary to take arbitrary linear combination (j){x) = A(j)i{x) + 502 (a;), A,B real 
constants, and look for solutions of the Schrodinger equation whose Wronskian with (j){x) vanishes 
for X —>■ 0: 

iy(0,^/>) = 0(x)^'(a;) -0'(x)^/'(a;) 0, x 0. (8.10) 

Physically this means that the particle's behavior near the point of singularity is fixed by the 
choice of auxiliary function (j){x). When (j){x) and its derivative take finite values at x = 0, the 
constraint ( p.lO[ ) becomes equivalent to the well known condition ip'{0) = ctp{0), where c is a real 
constant. 

For h < derivatives of all wave functions are singular at x = 0, and one needs to use a 

limiting procedure. Let us take as auxiliary function the general zero energy eigenfunction 



0(x) = ai/r cos(crlna; + 6'), a = y — 1/4, 
where a and 9 are some real constants. Substituting it into the above relation one finds, 

Vatan(alna; + ^)-^ ^ 0, a; ^ 0. 



(ilnx \/x \fx 



The chosen auxiliary function is homogeneous under the scaling of coordinate by a specific constant 
q = e^l'^ , (piix) = —^/q(p{x). Therefore, from the scaled form of the condition ( p.lO|) 

(j){qx)ilj' (qx) — (f)'{qx)tlj{qx) oc (j){x) '^^^^'^^ _ ^(^g^-^ _^ 

dx dx 



52 



it is seen that if iIj{x) satisfies ( |8.10|) , the same holds for the scaled wave function iplqx), i.e. scaling 
by g is a physical symmetry of the problem. This fact guarantees that discrete eigenvalues appear 
in the form of a geometric series. Since the behavior of the general solution of the Schrodinger 
equation for small x is known, the Wronskian vanishes only if near zero ip{x) has the form of (f){x) 
taken with the same angle 6, which is a free parameter of the self-adjoint extension. The spectrum 
itself is found from the requirement of normalizability of wave functions with such a property. The 
described technique of fixing the boundary condition at zero seems to be valid in the arbitrary 
case. 

As it was remarked in 0, eigenfunctions of the Schrodinger operators with self-similar poten- 



tials resemble discrete wavelets - the functions ip{x) affine transformations of which, ijjj^n{x) = 
2^/'^ip{2^x + n), generate an orthonormal basis of Hilbert space L^(R), {'4'j,n\''Pi,m) = ^ji^nm- In the 
above g-Floquet theory, we have a partial realization of this construction because the scaling of 
the coordinate of the wave functions by q creates an orthogonal function. The connection with 
wavelets originates, of course, from the group of affine transformations, because wavelets can be 



interpreted as coherent states associated with unitary representations of this group |^ according 
to one of the definitions mentioned in the introduction. An interesting fact that differential-delay 
equations similar to those determining coherent states of self-similar potentials may have solutions 



with finite support, or atomic solutions |39| indicates also a hidden relation with wavelets. 



9. The factorization method and "quantum Galois the- 
ory" 

Let us discuss the problem of integrability of a given equation. If this equation is algebraic, 
P{x) = Z)j=o'^i^"' — O5 with coefficients aj from some fixed number field k, then its solvability in 
terms of the radicals is determined by the solvability of its Galois group (the group of permutations 
of roots of P{x)) [0. One can say that algebraic equation is exactly solvable when the latter 
situation takes place, but this criterion is sensitive to the choice of k - there exists a universal 
field of complex numbers for which any algebraic equation is exactly solvable. 

Integrability, or exact solvability of differential equations is a looser notion. The weakest def- 
inition requires the existence of a solution of an initial value problem analytical in a bounded 
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domain. The class of equations integrable in this sense is too large - the field of functions deter- 
mined by formal power series can be thought as a universal one since solutions of any differential 
equation with coefficients from this field are given by functions of the same type. Another ex- 
treme definition consists in demanding for solutions to belong to differential fields of elementary 
or classical special functions, i.e. the fields built from the simplest functions and their derivatives. 
The question of integrability acquires precise meaning when the functions entering the equations 
and those allowed in the solutions are completely specified. In short, the problem is to understand 
transcendence of solutions over the given field of functions. 

In fact there exists a natural analog of the Galois theory for fields with operators |95| , |96| , which 



can be considered physically as a "quantum" (first-quantized) Galois theory. The main notions of 
it are the Picard-Vessiot (P-V) extension and the differential Galois group. Let u{x) in y"{x) = 
u{x)y{x) belong to the differential field k =C(a;) (field of rational functions). Roughly, the P-V 
extension M = k{yi, y2) is built from rational functions of independent solutions of this equation 
and their derivatives with coefficients from k. The group of automorphisms of k{yi,y2) which 
commute with the derivative and keep k fixed is called the Galois group Gal{M/k). For a general 
nth order linear ODE it is isomorphic to the group of invertable matrices GL{n,C). Roughly, the 
equation is solvable in terms of Liouvillian functions (exponentials, integrals, algebraics of the 
functions from k) when Gal{M/k) is solvable (or more simply when all elements of it can be made 
triangular). This theory is not widely known although it forms the basis of the intensively used 
computer programs of symbolic integration of differential equations [p^ . 

Let D be an indeterminate over some number field, then permutations of roots of a polynomial 
of D preserve this polynomial {D — xi) . . . {D — x„) = ^ {D — SjJ . . .{D — Xi^) = 0. In the 
quantum picture, when D is an operator, e.g. d/dx, there are two analogs of these "classical" 
permutations - the automorphisms figuring in the definition of the Galois group, which preserve the 
differential equation, and the permutations fi{x) ^ fj{x) in {D — fi{x)) . . . {D — fn{x))y{x) = 0. 
The latter do not, in general, preserve the initial equation, but cyclic permutations of such type 
play a key role in the factorization method H, |9^, which thus should be considered as one of the 



ingredients of "quantum" Galois theory. 

Consider the free Schrodinger equation with zero spectral parameter 



y" = {D + l/x) {D - l/x) y = (9.1) 
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over the field of complex numbers C. The P-V extension coincides with C{x) and the Galois 
group consists of nontrivial triangular matrices. Permutation of operator factors in ( |9.1| ) gives the 
equation 

{D - 1/x) {D + 1/x) y = y"- {2/x')y = 0, (9.2) 

which is different from the initial one, e.g. the potential lies now in C{x). Its Galois group is 
trivial (consisting of the unit matrix), because the general solution y = a/x + bx^ belongs to C(x). 
In general, the Darboux transformation changes the field of coefficients of ODE and its Galois 
group. However, an important feature is preserved in the chosen example; namely, solutions of 
both equations belong to one and the same field. If one factorizes the second equation using its 
general solution and permutes the operator factors, the solutions of the third equation, 

(O' + ^ff^) M.) ^ 0. (9.3) 

belong to C(x) again. This procedure of building rational potentials out of the zero one with the 
help of Darboux transformations with zero eigenvalue level has been considered in |10CI| , |101 



It can be shown that solutions of all equations built in this way lie in the P-V extension of ( p.l| ), 
i.e. the corresponding Galois groups are always trivial (an analogous situation takes place for 
repeated Darboux transformations with non-zero eigenvalue level). A similar iso-Galois picture 
prevails, e.g. for the Bessel functions. It would be interesting to investigate from this point of 
view the finite-gap potentials for which one has a simple factorization of the Hamiltonian. 

It is not difficult to describe a differential field which is preserved by any Darboux transfor- 
mation. For this, consider spectral problems generated by two Hamiltonians neighboring in the 
factorization chain 

{D ± f{x)){D T /(x))yi,2(x) = Xyi,2{x). (9.4) 
Let y{x, A) be a particular solution of the first equation. Then 

ai + hj ^^(^T^j (9-5) 

is a general solution. Consequently, y2{x, X) = {D — f{x))yi{x, A) gives a general solution of the 
second equation for A 7^ 0. For A = one has 

2/2(x, 0) = (^a2 + h f y\x\ 0)rfx') . (9.6) 



55 



This function belongs to the P-V extension of the first equation when the derivative dy{x, X)/dX 
at A = does so. Indeed, one can verify by direct differentiation that 

ry\x', \)dx' = y\x, + const. (9.7) 

J dxdX 

This relation shows that Darboux transformations preserve the differential field generated by 
solutions of a spectral equation and their derivatives over A for arbitrary values of A. This is a 
very wide field, but it exists and may be used in a search for such potentials, m(x), and values of A 
that ?/(x. A) are Liouvillian over the field from which u{x) is taken. It seems that the differential 
Galois theory, in conjunction with the equivalence up to a change of variables, can be used for 
a constructive definition of the notion of self-similarity. Indeed, for u{x) = and fixed A the 
P-V extension contains only exp{±n^/Xx)^ n gN, and this hints immediately at the form of 
the spectrum, A„ oc Ara^, for which Darboux transformations preserve this field. Analogously, 
defining relations of the self-similar potentials resemble the requirement of preservation of the 
P-V extension after N Darboux transformations. 

In a sense, the power of the factorization method (or of Darboux transformations) stems from 
its relation to the Galois theory. This explains also a large number of its applications. For 



instance, we mention its effectiveness in numerical calculations ||102||, in approximation theory 



103|| , in search of bispectral equations [|104|| , in the theory of orthogonal polynomials |59| , 



its connection with symplectic structures ||105|| , etc. In particular, Darboux transformations with 



zero eigenvalue appeared to play an important role in the construction of equations satisfying the 
Huygens principle ||100| , |1Q6|| . 



In conclusion, let us summarize the criteria which allow to call us potentials "exactly solvable" . 
First, such potentials should belong to some sufficiently simple differential field k, e.g. to C(x) but 
not to the field of formal power series. This requirement assumes that a full analytical structure 
of potentials is known. Then one can say that the Schrodinger equation is exactly solvable if 
the solutions satisfying the taken boundary conditions (i.e. when A belongs to the spectrum) 
are given by Liouvillian functions over k. This definition is tied to the differential Galois theory 
where one has a simple test of such solvability. Note that in this case a change of boundary 
conditions may change the "solvability" of the equation. A weaker definition of exact solvability 



refers to the availability of global structure of solutions [|107|| , i.e. it demands knowledge of various 



asymptotics of solutions, which often may be found even if the solution is not Liouvillian. The 
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latter requirement is natural for spectral problems, since in order to satisfy boundary conditions 
one should be able to connect solutions at various distances. In this respect, presently only classical 
special functions provide a completely satisfactory bank of information. 

The situation with self-similar potentials is instructive. It is difficult to specify the field of 
functions to which they belong (in particular, to find analytic properties of fj{x)), but once it is 
done the normalizable wave functions are given by Liouvillian functions. The discrete spectrum of 
these potentials corresponds thus to the exactly solvable problem in the differential Galois theory 
sense. Analysis of the continuous spectrum states, or of the solutions satisfying different boundary 
conditions, requires information which is not available at present, i.e. the rank of special functions 
for solutions of Schrodinger equations with self-similar potentials is not yet established. 

Coherent states provide a basis of the Hilbert space of states of quantum particle different 
from that determined by the Hamiltonian eigenstates. It would be interesting to understand their 
role from the differential Galois theory point of view. For the standard harmonic oscillator case 
the picture is simple. The second order differential (Schrodinger) equation is replaced by a first 
order one whose solutions provide an overcomplete set of Hilbert space vectors. Certainly this 
provides a simplification and "minimization" of the problem - the non-physical eigenfunctions 
of the Hamiltonian are removed in this procedure. Such minimization of the types of functions 
relevant to the given physical problem is characteristic for some coherent states from the functional- 
analytic point of view. This is evident for coherent states defined as orbits of states generated by 
physical symmetry groups. The ladder operator approach does not obey such a property in general 
because, starting from N = 2, coherent states of the self-similar potentials are determined by 
differential or differential-delay equations of higher order, which may contain nonphysical solutions 
as well. Moreover, for g 7^ 1 solutions of the latter equations contain some functional arbitrariness 
with respect to the Schrodinger equation. 



10. Conclusions 



The conclusions are short. The superpositions of coherent states in the abstract form ( p.l4|) , 
generalizing the even and odd states of |]r2|, are applicable to a very wide variety of systems. The 



parity coherent states (2^), ( 2.17 ) are less universal - their meaning and the way of derivation are 
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strongly tied to the presence of the parity symmetry. For any parity invariant potential admitting 
ladder operators A"!", A, one can make the canonical transformation A VA, A^ A'^V, where 
is a unitary operator linear in the parity operator (P]TH), which does not change the algebra 
satisfied by these operators and the Hamiltonian. This is the most simple physical observation 
of the present work. For the discussion of experimental implementations of superpositions of 
coherent states we refer to the papers |]ll], ^ |2^, and references therein. 



Coherent states of the free particle, constructed in section 4, differ qualitatively from those of 
the harmonic oscillator. First, they are built from the continuous spectrum states, which does 
not allow us to form superpositions similar to ( p.l4| ). Second, their analytical properties are quite 
unusual, so that at present it is not known even whether they are complete or not. It would be 
interesting to investigate the physical characteristics of these states. 

In many respects the results of this paper are not complete. This is caused by the complicated 
structure of self-similar potentials and their limiting cases, namely, by the appearance of Painleve 
transcendents and their g-analogs. However the author believes that this is a temporary situation 
and later many things will be simplified and, more importantly, the appropriate phenomenological 
applications of the coherent states for self-similar potentials considered in sections 5 and 6 will be 
found. Because of such expectations, the open problems indicated in this paper are interesting 
primarily from the physical point of view. Among the problems worth further investigation, we 
mention the detailed analysis of various limiting cases of the type discussed in section 7, the q- 
Floquet theory, connections with wavelets, etc. On the mathematical side, it is necessary to find 
general structure of g-special functions behind the self-similar potentials and to understand in a 
more appropriate way the differential Galois theory origin of the notion of self-similarity. 
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